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INTRODUCTION 


The  text  of  this  technical  report  is  very  brief  because  there 
is  no  better  way  to  describe  our  progress  in  two  years  than  by  the 
papers  we  have  written.  These  are  appended  to  this  report, 
separated  by  a  blank  colored  page.  Since  two  years  is  too  short  a 
tine  to  do  experiments,  analyze  the  data,  write  papers,  and  wait 
for  the  publication  process,  most  of  these  are  in  preprint  form. 
The  main  body  of  this  report  therefore  begins  with  several 
abstracts,  then  there  are  two  published  papers,  and  then  5 
preprints.  All  of  these  papers  have  been  submitted,  and  there  are 
still  others  in  progress  that  will  also  acknowledge  AFOSR  support. 


LIST  OF  ATTACHMENTS 


1.  Collection  of  abstracts 

2.  Copies  of  published  papers  on  Seam  Profile  Flattening 
(Opt.  Lett.)  and  on  Transient  Laser  Cooling  (PRL) . 

3.  Preprint  of  Quantum  Calculations  paper  (accepted,  PRA,RC) 

4.  Preprint  of  Simplest  Laser  Cooling  System  paper  (JOSA) 

5.  Preprint  of  Bichromatic  Cooling  paper  (PRL) 

6.  Preprint  of  Quantum  Effects  in  Helium  Cooling  paper  (PRL) 

7.  Preprint  of  Isotropic  Light  Cooling  paper  (PRA) 
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Gupta,  Padua,  Xie,  Batelaan,  and  Metcalf:  Longitudinal  Fields. 


Longitudinal  Magnetic  Field  Effects  in  Sub-Doppler  Cooling 
R.  Gupta,  S.  Padua,  C.  Xie,  H.  Batelaan,  and  H.  Metcalf 
Physics  Dept.,  S.U.N.Y.,  Stony  Brook,  NY,  11790 
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The  experimental  discovery  of  sub-Doppler  laser  cooling  was 
soon  followed  by  theoretical  descriptions  that  depended  on  optical 
pumping  processes  among  the  sublevels  of  multiply  degenerate 
atomic  ground  states.  These  theories  were  based  on  polarization 
gradients  in  count erpropagating  light  beams.  Later  it  was  found 
that  sub-Doppler  cooling  could  also  be  achieved  in  the  absence  of 
polarization  gradients  when  a  weak  transverse  magnetic  field  was 
used  for  mixing  eigenstates  of  the  light  shift  operator. 

By  contrast,  Zeeman  shifts  from  magnetic  fields  along  the 
quantization  axis  can  add  or  subtract  with  the  light  shifts  of 
those  states.  B{  thus  provides  an  additional  degree  of  freedom 
for  studying  magnetic  effects  in  laser  cooling  since  the  set  of 
sinusoidal  potentials  that  derives  from  the  spatially  periodic 
light  shifts  in  a  standing  wave  can  be  manipulated  with  the  field. 
For  example,  with  circularly  polarized  light  a  finite  Bf  can 
switch  the  normal  heating  to  cooling,  and  vice  versa.  We  have 
done  such  experiments  in  the  F  ■  3  <•  4  transition  of  85Rb  at  x  ■ 
780  nm  with  diode  laser  light,  and  have  shown  that  sub-Doppler 
heating  and  cooling  can  be  interchanged  with  opposite  longitudinal 
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field  directions,  in  good  agreement  with  this  model.  Such  experi¬ 
ments  may  also  allow  us  to  characterize  precisely  the  magnetic 
fields  seen  by  the  atoms  along  their  trajectories  in  the  beam. 

The  simplest  angular  momentum  schemes  that  permit  sub-Doppler 
cooling  have  4  levels,  and  one  example  is  an  F  >  1  *  o  transition 
in  either  a  linearly  or  circularly  polarized  standing  wave  with  a 
weak  transverse  magnetic  field.  The  28%  abundant  87Rb  isotope  has 
an  F  »  l  *  0  transition  that  is  ideal  for  studying  this  simplest 
possible  case.  Magnetically  induced  laser  cooling  can  reduce 
atomic  kinetic  energies  to  below  the  peaks  of  the  spatially  peri¬ 
odic  potential  seen  by  the  atoms  as  a  result  of  their  light 
shifts,  thus  channeling  them  between  the  planes  of  the  standing 
wave.  Furthermore,  their  deBroglie  wavelengths  become  comparable 
to  the  channel  size  so  that  their  motion  is  quantized  in  these 
ha If -wave length  sized  optical  traps. 

This  F  »  1  <*  0  transition  produces  only  one  light -shifted 
level,  and  thus  only  one  ground  state  sublevel  that  can  support 
bound  quantized  states.  An  rf -induced  transition  that  would 
change  the  populations  of  the  quantum  vibraticnal  states  of  this 
potential  (see  Fig.  1)  would  be  particularly  easy  to  interpret, 
and  would  be  subject  to  easily  understood  Zeeman  shifts  in  a 
longitudinal  field.  We  are  presently  studying  this  problem. 

Supported  by  NSF,  ONR,  AFOSR,  and  CAPES  (Brazil) 
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FIGURE  CAPTION 


Figure  1.  Energies  of  quantized  states  in  the  optical  light  shift 
potential  for  an  F  «  1  *  0  transition  in  circularly  polarized 
light  with  no  magnetic  field  (left)  and  positive  longitudinal 
magnetic  field  (right) . 


I 


Abstract  Submitted 
for  the  May  1993  Meeting  of  the 
Division  of  Atomic,  Molecular  and  Optical  Physics 


Physical  Review 
Analytic  Subject  Index 
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A  Simplest  Possible  Scheme  E&I  Sub-Doppler  laser 
Cooling.  R.  GUPTA,  S.  PADUA,  C.  X1E,  H.  BATELAAN,  and  H. 
METCALF.  SUNY  Stonv  Brook  -  The  simplest  angular  momentum 
schemes  that  permit  sub-Doppler  cooling  have  4  levels,  and  one 
example  is  an  F  «  1  *  0  transition  in  either  a  linearly  or  cir¬ 
cularly  polarized  standing  wave  with  a  weak  transverse  magnetic 
field.  The  28%  abundant  ,7Rb  isotope  has  an  F  **  1  -  0  transi¬ 
tion  that  is  ideal  for  studying  this  case.  Magnetically  induced 
laser  cooling  can  reduce  atomic  KE  to  below  the  peaks  of  the 
spatially  periodic  potential  seen  by  the  atoms  as  a  result  of 
their  bght  shifts,  thus  channeling  them  between  the  planes  of 
the  standing  wave.  Furthermore,  their  deBroglie  wavelengths 
become  comparable  to  the  channel  size  so  that  the  motion  is 
quantized  in  these  half-wavelength  sized  optical  traps.  This  F 
=  1*0  transition  produces  only  one  light-shifted  level,  and 
thus  only  one  ground  state  sublevel  that  can  support  bound  quan¬ 
tized  states.  An  rf-induced  transition  that  would  change  the 
populations  of  the  quantum  vibrational  states  of  this  potential 
would  be  particularly  easy  to  interpret,  and  would  be  subject  to 
easily  understood  Zeeman  shifts  in  a  longitudinal  field.  We  are 
presently  studying  this  problem. 
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Raman  Processes  in  Sub-Doppler  Laser  Cooling,  c.  XTc,  s. 
PADUA,  R.  GUPTA,  H.  BATELAAN,  and  H.  METCALF,  SUNY 
Stony  Brook  -  We  have  observed  a  new  type  of  sub-Doppler  laser 
cooling  with  neither  polarization  gradients  nor  magnetic  fields. 
All  known  sub-Doppler  laser  cooling  methods  require  one  of  these 
two  to  redistribute  the  atomic  population  among  the  energy  le¬ 
vels.  This  allows  the  irreversible  processes  of  optical  pumping 
(OP)  to  manipulate  the  conservative  force  of  light  shifts  (e.g. 
resonant  exchange  between  the  two  beams  comprising  a  standing 
wave)  [1]  into  a  damping  force.  Our  scheme  employs  stimulated 
Raman  transitions  between  ground  state  hyperfine  levels  of  Rb 
with  a  bichromatic  standing  wave,  accompanied  by  a  momentum 
change  of  Ap  »  2hk.  OP,  with  Ap  «  1  hk,  creates  a  population 
imbalance  allowing  repetition  of  the  Raman  process.  When  the 
detuning  of  the  light  from  Raman  resonance  5  is  properly  chosen, 
the  net  momentum  exchange  damps  the  atomic  motion.  Changing  6 
by  only  jr/6  (1  MHz)  replaces  cooling  by  heating.  This  process 
separates  the  velocity  selective  Raman  resonance  and  OP  proces¬ 
ses  differently  from  others,  and  helps  clarify  the  role  of  these 
different  mechanisms. 

1.  G.  Nienhuis  et  al.,  Phys.  Rev.  A44,  462  (1991). 
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Transient  Effects  in  Lassi  Cooling.  S.  PADUA,  C.  XIE,  R. 
GUPTA,  H.  BATELAAN,  T.  BERGEMAN,  and  H.  METCALF, 
SUNY  Stonv  Brook  -  Transient  laser  cooling  (TLC)  can  produce 
cooling  and  heating,  but  often  with  the  opposite  detuning  from 
that  found  in  steady  state.  In  TLC  the  time  scale  is  set  by  the 
optical  pumping  (OP)  rate  to  a  state  not  coupled  by  the  laser 
field.  The  combination  cf  such  OP  processes  and  the  conservative 
light  shift  potential  UpSinHz  leads  to  TLC.  The  average  PE  of 
atoms  entering  a  standing  wave  is  Uq/2.  They  experience  the  op¬ 
tical  force  until  undergoing  OP  to  an  uncoupled  state,  which  is 
more  likely  to  happen  at  high  light  intensity,  near  an  antinode. 
For  a  >  0  this  means  higher  PE  and  thus  lower  KE,  and  conversely 
for  8  <  0.  In  TLC  there  is  no  final  “temperature"  resulting  from 
competition  between  a  damping  force  and  diffusive  heating.  In¬ 
stead  the  changes  in  KE  are  bounded  by  U0  so  that  the  signal 
widths  decrease  with  intensity.  This  can  result  in  sub-Doppler 
widths.  We  have  made  two  independent  theoretical  studies  of  these 
experiments.  In  a  semiclassical  calculation  we  evolve  the  motion 
for  r  calculated  OP  time  and  calculate  the  velocity  distribution. 
We  have  also  performed  fully  quantum  mechanical  calculations  of 
the  motion  of  atoms  in  the  standing  wave  whose  basis  set  consists 
of  product  states  of  internal  and  external  atomic  coordinates. 

Supported  by  NSF,  OMR,  AFOSR,  and  CAPES  (Brazil) 
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Longitudinal  Magnetic  Field  Effects  ip  SubDnppler 
Cooling.  R.  GUPTA,  C.  XIE,  S.  PADUA,  H.  BATELAAN,  and  H. 
METCALF.  SUNY  Stony  Brook  --  Sub-Doppler  laser  cooling  car.  be 
achieved  in  the  absence  of  polarization  gradients  when  a  weak 
transverse  magnetic  field  is  used  for  mixing  eigenstates  of  the 
light  shift  operator.  By  contrast,  Zeeman  shifts  from  magnetic 
fields  B*  along  the  quantization  axis  can  add  or  subtract  with 
the  light  shifts  of  those  states.  B$  thus  provides  an  addi¬ 
tional  degree  of  freedom  for  studying  magnetic  effects  in  laser 
cooling  since  the  set  of  sinusoidal  potentials  that  derives 
from  the  spatially  periodic  light  shifts  in  a  standing  wave  can 
be  manipulated  with  the  field.  For  example,  with  circularly 
polarized  light,  a  finite  B{  can  switch  the  normal  heating  to 
cooling,  and  vice  versa.  We  have  done  such  experiments  in  the 
F  =  3  *  4  transition  of  °5Rb  at  X  *  780  nm  with  diode  laser 
light,  and  have  shown  that  sub-Doppler  heating  and  cooling  can 
be  interchanged  with  opposite  longitudinal  field  directions,  in 
good  agreement  with  this  model.  One  application  of  such 
experiments  is  to  characterize  precisely  the  magnetic  fields 
seen  by  the  atoms  along  their  trajectories  in  the  beam. 
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SufcDflfiact  User  Coding  of  Mstastafrle  Hs*  M.  widmer, 

M.J.  BELLANCA,  E.  VREDENBREGT,  T.  BERGEMAN,  and 

H.  METCALF,  SUNY  Stonv  Brook  -  We  have  observed  magnetic-  1 

ally  induced  laser  cooling  (M1LC)  on  the  J  «  t  *  2  component  of 

the  23S  •  23P  transition  at  x  *  1.083  im:  at  metastable  helium 

(He*).  At  low  magnetic  field  the  measured  widths  of  our  atomic 

velocity  distributions  are  as  narrow  as  2hk>  corresponding  to  a 

temperature  »  3  pK  (1-D  Doppler  limit  a  25  X).  At  higher  fields 

(-150  mG)  the  velocity  distribution  split,  into  two  peaks  of  1 

sub-Doppler  width  whose  separation  is  prop:  tional  to  the  Zeeman 
splitting.  Our  measurements  are  cons** '  with  model  calcula¬ 
tions  based  on  the  time  evolution  of  jantum  density  matrices. 

Our  experiments  explore  the  transition  regime  between  the  quan¬ 
tum  mechanical  and  semiclassica!  m  '  :els  of  momentum  exchange  in 

atom-light  interactions  because  oi  the  large  ratio  (0.4)  of  the  1 

recoil  velocity  tik/M  to  the  rms  of  the  Doppler-limited  velocity 

distribution  ih»/3M  for  this  transition.  We  use  a  supersonic 
LN,  cooled  nozzle  source  excited  by  a  de  discharge  for  He*,  and 
we  excite  the  23P  state  with  a  home-built  diode-laser-pumped  LNA 

laser.  I 

Supported  by  NSF,  ONR.  AFOSR,  and  Dept.  Educ. 


(  )  Prefer  Poster  Session 
(  )  Prefer  Standard  Session 
(  )  No  Preference 


Submitted  by  Harold  Metcalf 


js  jr .  *s  >; 

Stpaui 

‘  / 


Harold  Metcalf 
Physics  Depan  men? 

Sate  University 
Stony  Brook.  NY  i\m 


» 


* 


•  •  •  •  •  i 


* 


Abstract  Submitted 
for  the  May  1993  Meeting  of  the 
Division  of  Atomic,  Molecular  and  Optical  Physics 


» 


» 


Physical  Review 
Analytic  Subject  Index 
Number:  32.80.Pj 


Bulletin  Subject  Heading 
in  which  paper  should  be  placed: 
Laser  Cooling 


Hr) 


Quantum  Qalguiatims  £21  One-Dimensional  Cooling  Qf 
Helium.  E.  VREDENBREGT,  M.  DOERY,  T.  BERGEMAN  and 
H.  METCALF,  SUNY  Stony  Brook  -  We  report  theoretical  velocity 
distributions  for  sub-Doppler  laser  cooling  of  metastable 
He*(23S),  calculated  with  the  Density  Matrix  l)  and  Monte  Carlo 
Wavefunction  2>3>  approaches.  For  low-field  (B  =  50  mG) 
magnetic-field  induced  laser  cooling  on  the  23S  -»  (23P,  J  =  2) 
transition  (x  =  1083  nm),  we  get  a  narrow,  sub-Doppler 
structure,  consisting  of  three,  -1  photon  recoil  wide  peaks, 
spaced  -1  recoil  apart.  With  increasing  field,  this  three-peak 
structure  develops  into  two  velocity-selective  resonance  (VSR) 
peaks,  each  -2  recoils  wide.  For  the  23$  -*  (33P,  J  =  2) 
transition  (A  =  389  nm),  VSR  peaks  are  predicted  to  appear  at 
low  field  without  the  third,  central  peak,  which  only  develops 
at  higher  field  (B  =  200  mG).  Additional  computations  deal  with 
polarization-gradient  cooling.  In  genera],  we  find  that  for  one¬ 
dimensional  cooling  calculations,  the  Density  Matrix  method  is 
more  efficient  than  the  Monte  Carlo  Wavefunction  approach. 
Experiments  are  currently  under  way  to  test  the  results. 
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DARK  STATES  AND  LASER  COOLING 
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L  Introduction.  We  present  the  results  of  a  series  of  one-dimensional  laser  cooling 
experiments  where  atomic  dark  states  play  an  important  role.  In  some  cases  the  atoms  are 
removed  from  the  dark  states,  and  this  process  is  vital  to  the  cooling  mechanism.  In  other 
cases  the  dark  states  are  velocity  dependent,  and  thus  atoms  in  such  states  arc  "cc’d"  by 
virtue  of  their  narrow  velocity  distribution.  All  the  experiments  are  done  in  a  one-dimen¬ 
sional  standing  wave  at  sufficiently  low  intensity  and  large  detuning  that  the  excited 
state  population  is  small.  The  laser  beams  are  crossed  by  an  atomic  beam  that  is  cooled  in 
its  transverse  direction  by  its  interaction  with  the  light  field. 

We  define  a  generic  dark  state  as  a  stationary  state  of  the  Hamiltonian  tf  that  can  not 
be  excited  by  the  light  field.  Jf  must  include  the  atom-light  interaction  as  well  as  the 
kinetic  energy  of  Die  atomic  motion.  In  principle,  atoms  in  a  dark  state  could  be  detected 

by  magnetic  deflection  or  by  excitation  from  a  probe  laser  beam  not  included  in  N. 

* 

ILL  Transient  LiSEl  Cooling.  In  the  first  of  the  four  studies  reported  here,  we 
consider  cooling  in  the  non-steady  state  case.  The  optical  Bloch  equations  cannot  be 
solved  by  setting  dp/dt  «  0  as  is  usually  done,  because  now  atoms  are  quickly  pumped  into 
uncoupled  ground  states  (e.g.  alkali  hfs)  at  a  rate  rp.  Even  in  the  short  optical  pumping 
time  1/yp,  atoms  are  subject  to  the  conservative  force  from  the  light  shift  potential 
(calculated  by  adiabatically  eliminating  the  excited  state).  Its  sinusoidal  spatial  depen¬ 
dence,  together  with  optical  pumping  to  the  dark  state,  can  be  exploited  to  reduce  the 
average  KE  of  the  atoms. 

We  consider  the  case  of  a  blue-tuned  User  field  whose  light  shift  is  positive.  Atoms 
that  enter  in  the  antinodes  ire  quickly  pumped  to  the  dark  state  by  the  high  intensity 
light  before  there  is  enough  time  for  their  KE  to  change  significantly.  By  contrast,  atoms 
that  enter  in  the  nodes  must  travel  to  a  higher  intensity  region  before  they  can  be  pumped, 
and  in  doing  so  oust  work  against  the  positive  light  shift  and  therefore  be  slowed.  Thus 
on  the  avenge  atoms  are  slowed,  and  cooled.  For  red  tuned  light  the  atoms  are  heated. 
Once  pumped  to  the  dark  state,  the  cooled  atoms  are  unperturbed  by  cither  the  conservative 
force  or  momentum  diffusion  because  they  no  longer  interact  with  the  light. 

ILL  .Vcfaily  Selective  fitthcsal  EflguUliflfl  logging*  (VSCPT)  In  some  arrangements 
of  polarization  and  atomic  angular  momentum,  atoms  irradiated  by  light  beams  of  different 
polarizations  may  be  optically  pumped  to  a  coherent  dark  state.  This  occurs  if  the  light 
beams  are  resonant  with  a  Raman  transition  to  provide  coherence.  For  a  degenerate  pair  of 
ground  state  sublevels,  the  two  light  beams  must  have  the  same  frequency  in  the  atomic 
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frame.  When  they  also  have  the  same  frequency  in  the  lab  frame,  then  only  atoms  with  <  v  >  - 
0  in  a  coherent  superposition  states  of  v  =  ±hk  can  populate  the  dark  state.  Moving  atoms 
will  see  opposite  Doppler  shifts  of  the  two  laser  beams  so  the  Raman  resonance  condition  is 
not  fulfilled  unless  the  atomic  energies  (including  KE)  are  appropriately  split.  For  cases 
studied  up  to  now,  the  dark  state  can  only  be  reached  by  a  random  walk  in  momentum  space 
(associated  with  spontaneous  emission)  with  steps  of  size  hk.  Thus  the  rate  of  populating 
die  dark  state  is  quite  small. 

Recently  two  groups  have  suggested  that  the  polarization  scheme  lin  u  lin  at  45°  can 
both  produce  dark  states  and  provide  polarization  gradient  cooling  into  them  We  have 
observed  the  dramatically  increased  rate  of  populating  these  states  in  the  J  =  1  ■*  1  tran¬ 
sition  of  metastable  He.  This  has  enormous  promise  for  production  of  diffraction  limited 
atomic  beams  (limited  only  by  their  deBroglie  wavelength).  We  are  now  studying  this  effect 
for  quantitative  comparison  with  theory. 

the  pol^i^on^dlenfand^  magnetic  mixing  schemes^of  SDLC require  multiple  ground  state 
sublevels  and  optical  pumping  among  them.  The  natural  question  to  ask  is  "What  is  the  sim¬ 
plest  possible  level  structure  that  permits  all  types  of  SDLC?"  For  realistic  situations, 
this  "hydrogen  atom  of  SDLC"  is  an  Fg  =  1  *  Fc  =  0  transition.  Such  a  system  always  has 
dsurk  states  in  1-D,  although  the  choice  of  basis  determines  their  composition  from  the  Mp 
states.  In  a  pure  standing  wave  there  are  always  two  dark  states,  but  in  a  polarization 
gradient  there  is  only  one.  Atoms  pumped  into  these  may  be  ^leased  firm  them  by  an 
appropriate  B  field. 

We  have  demonstrated  all  the  different  types  of  SDLC  in  this  system  using  ^Rb.  These 
include  lin  4  lin  and  <r+  >  <r*  polarization  gradients  in  a  B  field  to  empty  the  dark  states. 
The  B  field  must  not  be  parallel  to  any  polarization  axes,  and  must  be  strong  enough  to 
define  the  quantization  axis.  All  these  experiments  can  be  described  by  a  Sisyphus  picture 
of  cooling  to  v  =  0.  Furthermore,  we  have  observed  several  examples  of  SDLC  to  a  non-zero 
velocity  proportional  to  |B|  using  velocity  selective  Raman  resonances  (VSR).  All  our  F  = 
1  *  0  experiments  can  be  described  in  the  momentum  exchange  picture  of  VSR. 

P,4.  Two  Frequency  SDLC.  Selection  rules  for  the  F  =  1  *  0  excitation  do  not  allow 
decay  to  a  dark  hfs  state,  so  re-pumping  is  not  necessary,  but  this  is  not  so  for  all  exci¬ 
tations.  In  <”Rb  we  used  two  frequencies  of  light  generated  by  an  EOM  and  tuned  to  couple 
the  hfs  ground  states  via  Raman  resonance.  When  the  nodes  of  the  two  standing  waves  coin¬ 
cide,  there  is  only  Doppler  cooling.  However,  when  the  nodes  are  separated  by  A/4  there  is 
another  type  of  Sisyphus  cooling  with  both  lasers  tuned  blue  of  resonance,  or  heating  witn 
red  tuning.  We  have  measured  the  dependence  of  the  cooling  efficiency  on  this  spatial 
phase  shift  and  found  it  sinusoidal,  as  expected.  Furthermore,  we  also  observed  VSR  peaks 
that  shift  linearly  with  the  detuning  from  Raman  resonance  in  this  configuration. 

When  the  two  lasers  have  opposite  detuning,  one  red  and  one  blue,  a  deflection  force 
appears  that  does  not  change  sign  with  velocity.  We  measured  the  velocity  dependence  and 
capture  range  of  this  optical  rectification  force  and  found  it  to  be  consistent  with  a 
simple  model. 

III.  Conclusion.  In  sections  II.  1  and  II.2  we  described  experiments  where  atoms  are 
pumped  into  dark  sjates,  and  in  11,3  and  II.4  we  described  experiments  where  dark  states 
were  emptied  by  a  B  field  or  by  light  of  a  different  frequency.  In  each  example  atoms  were 
cooled  below  the  Doppler  temperature  (to  the  recoil  limit  in  VSCPT).  In  all  cases  there  is 
a  conservative  force  from  the  spatially  varying  light  shift  potential,  and  this  is  mani¬ 
pulated  into  a  dissipative  force  by  the  irreversible  process  of  spontaneous  decay. 

Supported  by  N.S.F.,  O.N.R.,  and  A.F.O.S.R. 
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We  describe  a  novel  technique  to  make  a  Gaussian  laser  beam  profile  spatially  flat.  We  exploit  the  angular 
dependence  of  the  transmission  of  an  6 talon  to  tailor  the  spatial  profile  to  the  desired  form.  A  simple  analysis 
shows  why  our  method  works  so  well  and  how  an  £  talon  could  be  tuned  to  give  the  optimum  results  at  all 
wavelengths.  This  technique  has  enormous  advantages  over  other  methods. 


The  output  beams  of  most  lasers  are  not  spatially 
uniform  because  the  resonant  modes  of  their  op¬ 
tical  cavities  have  Gaussian  profiles.  The  output 
beams  of  single-mode  lasers  are  usually  in  the  fun¬ 
damental  TEMoo  mode  that  has  a  Gaussian  spatial 
profile.  Since  a  number  of  experiments  require  a 
flat  spatial  beam  profile,  it  is  important  to  be  able 
to  transform  to  this  kind  of  laser  beam,  and  several 
techniques  for  converting  these  Gaussian  beams  to 
others  that  have  flatter  profiles  have  been  developed. 
We  demonstrate  a  new  but  relatively  simple  approach 
that  has  the  advantages  of  being  highly  efficient, 
quite  flexible,  and  low  cost. 

There  are  at  least  two  obvious  ways  to  perform  this 
task.  One  is  to  block  all  but  the  central  region  of 
the  laser  beam  with  a  small  aperture.  The  resulting 
beam  intensity  has  a  central  maximum  and  decreases 
toward  the  edge,  but  the  decrease  is  not  significant  if 
the  aperture  size  is  small  compared  with  the  Gauss¬ 
ian  beam  half-width.  The  advantages  of  this  method 
ate  its  low  cost  and  ease  of  implementation.  The 
principal  disadvantage  of  this  method  is  the  loss 
of  light.  In  order  to  have  the  intensity  difference 
between  the  center  and  the  edge  be  less  than  5^, 
the  aperture  can  only  transmit  less  than  5%  of  the 
total  beam  power.  Another  disadvantage  arises  from 
diffraction  by  the  aperture’s  edges;  the  experimental 
region  must  not  be  separated  from  the  aperture  by 
too  many  times  its  diameter,  or  imaging  optics  must 
be  employed. 

Another  way  is  to  use  a  Alter  whose  transmis¬ 
sion  is  not  spatially  uniform.  Consider  a  nonuniform 
Alter  whose  transmission  7V)  increases  from  some 
value  0  <  1  at  the  center  (r  -  0)  to  unity  at  some 
distance  r,  from  -he  center.  We  define  G/(r)  as  a 
Gaussian  profile  of  chare-teristic  width  tty  and  let  the 
transmission  vary  rs  l\r)  -  0/Gfir)  for  r  <  r„  where 
G/(rt)  -  0.  and  Ifr)  -  1  for  r  >  r„  The  profile  of  a 
Gaussian  laser  beam  of  waist  mi  -  tty  transmitted  by 
such  a  filter  is  uniform  over  the  region  from  r  -  0 
to  r,.  For  0  -  1,  rt  -  0,  and  the  power  of  such  a 
beam  is  zero,  but  for  0  -  1/e2,  this  area  extends 
until  rf  -  wi.  The  intensity  in  the  uniform  part  of 
the  beam  profile  is  0  times  smaller  than  the  peak 
intensity  of  the  original  Gavaian  beam,  and  the 
total  power  transmitted  in  the  fiat  profile  region  is 
-0in(0)  times  the  laser  beam  power. 


Such  a  filter  overcomes  the  great  loss  of  laser 
power  that  occurs  when  one  uses  an  aperture  as 
described  above,  but  there  are  at  least  three  severe 
disadvantages  to  this  method  of  beam  profile  flat¬ 
tening.  First,  although  such  filters  are  commercially 
available,  they  are  expensive  and  usually  custom 
made.  Second,  it  is  necessary  to  choose  a  filter  for 
a  particular  value  of  0  and  tty  and  then  match  the 
laser  beam  waist  to  the  filter.  Any  flexibility  would 
demand  an  assortment  of  such  filters.  Third,  mount¬ 
ing  and  alignment  of  the  filter  is  quite  critical;  it  must 
be  centered  precisely  on  the  laser  beam. 

We  have  developed  a  beam  profile  flattening  scheme 
that  has  many  of  the  advantages  of  both  of  these 
methods  without  their  disadvantages.  We  exploit 
the  angular  dependence  of  the  transmission  charac¬ 
teristic  of  a  thin  Fabry- Perot  6talon  to  tailor  the 
spatial  profile  of  the  beam  as  desired.  TVansmission 
of  a  nonparallel  beam  through  an  6talon  gives  the 
familiar  bull’s-eye  pattern  similar  tQ  Newton’s  rings. 
With  the  proper  choice  of  various  parameters,  the 
central  region  can  have  transmission  0  <  1  at  r  =  0, 
and  its  increase  toward  the  first  bright  ring  at  r  -  re 
can  vary  nearly  as  0 /G/<r). 

For  a  parallel  beam  of  light  whose  wavelength  is 
A  in  air  incident  at  angle  6 ,  the  transmission  in  air 
of  a  solid  6talon  of  thickness  d,  reflectivity  R,  and 
refractive  index  n'  is1 


7X0') 


1 

1  +  F  sins(o/2)  ’ 


(1) 


where  a  ■  4 nn'd  cos  O'/A  ■  a  cos  O’,  F  ■  4/f/(l  - 
R)3,  and  n'  sin  0'  ■*  sin  0.  This  function  is  plotted 
in  Fig.  1.  At  the  center,  where  0'  -  0  -  0,  the 
transmission  0  is  given  by  1/0  -  1  +  F  sinJ(a/2), 
and  0  <  1  as  expected. 

For  a  tunable  laser  in  an  application  where  the 
laser  frequency  is  not  crucial,  the  frequency  can  be 
adjusted  for  selection  of  any  chosen  value  of  0.  For 
an  4talon  ~T  mm  thick,  a  -  inn'd/X  -  10\  so  a 
shift  of  only  a  few  gigahertz  can  vary  (1  -  0)  by  a 
factor  of  3  or  more.  In  a  spectroscopy  experiment 
where  the  laser  frequency  is  constrained  by  atomic 
resonances,  the  value  of  0  can  instead  be  adjusted  by 
varying  d  or  rt'.  It  ia  particularly  convenient  to  use 
an  air-spaced  £  talon  that  can  be  tuned  by  varying  d 
with  piezoelectric  elements.  On  the  other  hand,  a 
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Fig.  1.  Angular  dependence  of  the  transmission  of  a 
Fabry -Perot  ttalon. 

solid  dtalon  can  be  tuned  by  varying  its  temperature, 
which  changes  o  through  both  n'  and  dP  For  an 
dtalon  of  modest  reflectivity  (F  -  3),  0  can  be  varied 
by  0.5  by  a  change  of  just  a  few  degrees  Celsius. 

lb  calculate  the  effect  of  the  filter  function  7X0') 
on  a  diverging  or  converging  Gaussian  laser  beam, 
we  note  that  because  a  is  so  large,  the  lowest-order 
region  requires  only  a  small  divergence  angle.  Thus 
we  can  restrict  our  attention  to  0'  «  1,  expand  cos  0', 
and  keep  only  the  first  two  terms.  We  write 

2  sina(a/2)  »  1  -  cos(a  cos  0') 

-  1  -  cos[a(l  -  0‘2/2  +  ...)] 

-  1  -  cos(a)cos(a0'2/2) 

-  sin(a)sin(a0'2/2)  +  . . . .  (2) 

Even  though  a  ~  104,  a0'2  varies  between  0  and  2 it, 
so  the  approximation  a0'4/ 24  «  1  in  Eq.  (2)  is  very 
well  satisfied. 

lb  see  that  such  a  filter  function  will  efficiently  flat¬ 
ten  the  profile  of  a  Gaussian  beam  we  note  that  for  a 
diverging  or  converging  laser  beam,  0  is  proportional 
to  the  distance  r  from  the  center  in  a  plane  perpendic¬ 
ular  to  the  beam  axis.  We  write  r  ■  0L,  where  L  » 
d  is  the  distance  from  the  beam  waist  to  the  dtalon, 
and  0  -  n'01  because  0  «  1.  We  write  the  beam 
profile  as  G,(r)  -  Gt(n'L6‘)  ■  exp[-2(n'L0'/w;)2],  set 
T(0')G|(0')  ■  0t  and  expand  both  T(0')  and  G;(0') 
about  0'  -  0.  We  find  that 

1  *-  2 (p0?  +  2(p0')4  +  ... 

-  0{1  +  (F/2)[l  -  cos(a)cos(a0'J/2) 

-  Bin(o)sin(a0'4/2)3} ,  (3) 

where  p  ■  n'L/wi  and  typically  p  »  1. 

Equating  the  0'-independent  terms  gives  0  as  dis¬ 
cussed  above,  and  choosing  a  particular  spacing  or 
temperature  sets  both  a  and  jS  for  a  given  dtalon 
(fixed  F).  Comparing  the  next-order  terms  in  0', 
we  find  that  the  condition  for  cancellation  of  the 
6*  term  is  0aF  sin(o)  -  8p\  and  it  can  be  satisfied 
by  varying  L  (and  thus  p)  with  a  small  telescope. 
Similarly,  the  condition  for  cancellation  of  the  6* 
term  is  fia*F  cos(a)  -  32p4.  This  last  condition  can 
be  met  by  selecting  F  and  modifying  the  italon  spac¬ 
ing  or  temperature  accordingly  to  maintain  the  de¬ 
sired  a.  The  two  conditions  together  require  that 
eoe(o)  -  F/(F  +  2).  Keeping  terms  in  0*  in  Eq.  (2) 
does  not  change  these  conclusions  but  only  causes 
small  changes  in  the  numerical  factors.  Figure  2 
ahows  plots  of  the  function  T{r)Gi(r)  for  a  range  of 


0’s.  Of  course,  the  choice  of  which  parameter  (T,  p, 
or  F)  to  cancel  which  term  in  0'  can  be  made  in  any 
order,  and  F  could  be  used,  for  example,  to  choose  0. 

This  system  has  several  obvious  advantages  over 
previously  used  beam-flattening  methods.  First,  it 
allows  maximum  flexibility  for  the  use  of  the  avail¬ 
able  power  and  still  achieves  extraordinary  unifor¬ 
mity  across  the  beam  profile,  as  shown  in  Fig.  2. 
Second,  it  is  easily  adjustable  for  beams  of  different 
sizes.  The  angular  divergence  is  simply  adjusted 
with  a  small  Galilean  telescope.  No  centering  of 
the  beam  or  the  btalon  is  needed.  Third,  its  cost 
is  much  lower  than  that  of  most  available  Gaussian 
filters,  and  suitable  dtalons  are  widely  available.  By 
contrast,  it  has  two  notable  disadvantages.  First,  0 
is  frequency  dependent,  and  in  a  scanned  frequency 
spectroscopy  experiment,  the  6talon  spacing  or  tem¬ 
perature  would  have  to  be  tuned  with  the  laser.  In 
addition,  there  might  be  some  problems  with  multi- 
mode  lasers,  but  since  typical  mode  separations  are 
a  few  hundred  megahertz,  the  effect  on  0  should 
be  small.  Second,  light  that  is  not  transmitted  by 
the  dtalon  is  reflected  back  to  the  laser,  where  such 
feedback  may  cause  stability  problems. 

In  our  experiments,  we  use  diode-laser  light  to 
excite  atoms  in  an  atomic  beam.  We  need  a  rect¬ 
angular  beam  of  light  that  is  uniform  only  along 
the  atomic  beam  direction.  Furthermore,  we  can 
afford  to  have  absolutely  no  reflection  back  to  the 
laser  because  diode  lasers  are  extremely  sensitive 
to  tiny  amounts  of  feedback.  Fortunately  the  6talon 
filter  described  above  offers  solutions  to  both  of  these 
problems  by  using  higher  orders  of  the  6talon  for 
beam  profile  flattening.  In  this  case,  the  beam  is 
diverging  in  only  one  dimension,  and  the  dtalon  is 
tilted  to  some  few-degree  angle  0„  with  respect  to 
the  laser  beam.  Thus  the  transmission  pattern  is 
a  series  of  unequally  spaced  parallel  stripes  instead 
of  rings,  and  we  work  in  one  of  the  outer  stripes. 
Of  course,  the  range  of  angles  0'  about  B0  is  still 
restricted  to  allow  sin(a/2)  to  range  over  less  than 
one  cycle. 

This  scheme  offers  two  additional  advantages  over 
the  zero-order  method  described  above.  First,  the 
choice  of  0O  offers  yet  another  free  parameter  to  help 
flatten  the  beam  profile.  For  example,  0„  can  be  var¬ 
ied  to  maintain  0  in  a  frequency-swept  spectroscopy 
experiment.  Second,  light  not  transmitted  by  the 
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Fig.  2.  Calculated  transmitted  profile  of  a  Gaussian 
beam  by  a  solid  ttalon  for  several  different  values  of  0. 
Here  0  is  varied  by  choosing  ft  to  be  10%  [curve  (a)},  30% 
[curve  <b)l,  and  60%  (curve  (c)). 
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Fig.  3.  Calculated  transmitted  profile  of  a  Gaussian 
beam  in  one  dimension  by  a  solid  6 talon  in  third  order 
(third  ring).  Curve  (a)  shows  the  transmission  function, 
curve  (b)  shows  the  Gaussian  profile,  and  curve  (c)  is  their 
product. 
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Fig.  4.  (a)  Measured  and  (b)  calculated  spatial  profiles 
of  a  diode-laser  beam  using  our  solid  4 talon. 

dtalon  is  reflected,  and  the  spatial  profile  of  the  light 
beam  reflected  out  at  angle  20e  has  a  strong  central 
intensity  maximum  that  can  be  easily  injected  into 
a  fiber  and/or  used  for  User  diagnostics  such  ss  sat¬ 
urated  absorption  for  laser  frequency  locking.  Note 
that  this  scheme  flattens  the  beam  profile  in  only  one 
dimension  instead  of  two,  but  that  is  precisely  what 
is  needed  for  exciting  atoms  in  a  beam  that  crosses 
the  laser  beam, 

The  few  center  stripes  of  the  transmission  pat¬ 
tern  are  rather  asymmetrical  about  their  minima 
(Fig.  1),  but  the  next  few  are  quite  satisfactory,  both 
because  their  asymmetry  is  relatively  small  and  be¬ 
cause  they  are  not  so  tiny  that  they  are  difficult  to 
use.  Expansion  about  9'  -  9,/n'  instead  of  0'  -  0  is 
straightforward,  but  in  that  case  fi  is  determined  by 
F  instead  of  being  s  free  parameter.  The  expansion 


is  not  discussed  here,  but  the  results  of  numerical 
calculations  are  presented  in  Fig.  3. 

We  have  measured  the  flattened  beam  profiles  pro¬ 
duced  this  way  for  several  different  laser  and  6talon 
parameters  and  present  the  results  in  Fig.  4.  In 
all  cases,  the  experimental  curves  in  Fig.  4(a)  are 
in  excellent  agreement  with  the  calculated  curves 
in  Fig.  4(b).  The  most  striking  aspect  of  Fig.  4  is 
not  the  good  agreement  between  the  measured  and 
calculated  beam  profiles  but  just  how  flat  the  profiles 
can  be  made  with  this  simple  device.  By  chang¬ 
ing  the  angle  of  the  4talon,  one  can  tailor  a  profile 
that  is  very  flat  from  one  end  to  uie  other.  The 
slight  asymmetry  introduced  by  the  asymmetry  of  the 
transmission  curve  can  be  partially  compensated  by  a 
slight  additional  tilt  and  can  be  minimized  by  using 
a  higher-order  transmission  curve.  We  found  that, 
over  a  wide  range  of  values  of  0,  typically  40%  of 
the  laser  beam  power  could  be  in  the  top-hat-shaped 
profile  this  way. 

All  the  measurements  were  done  with  a  Sharp 
LT021  5-mW  semiconductor  laser  operating  at  780 
nm  at  room  temperature.  The  solid  6talon  was  ss0.6 
mm  thick  and  coated  for  R  ss  30%  at  A  =  800  nm. 
The  spatial  profile  was  scanned  with  a  photodiode  by 
using  a  100-/im  slit  mounted  in  front  of  it. 

In  summary,  we  have  demonstrated  an  extremely 
flexible,  efficient,  and  low-cost  technique  to  make  a 
Gaussian  laser  beam  profile  spatially  fiat  by  exploit¬ 
ing  the  angular  dependence  of  the  transmission  of  an 
4talon.  The  method  can  be  easily  scaled  up  to  high- 
power  lasers  for  industrial  cutting  and  welding  and 
scaled  down  to  miniature  devices  for  printers,  com¬ 
munications,  and  medical  diagnostics  and  treatment. 
It  can  be  used  in  the  UV  for  curing,  photolithography, 
and  semiconductor  processing  and  in  the  IR  for  heat¬ 
ing.  It  is  applicable  to  any  Gaussian  beam  but  may 
be  especially  useful  with  diode  lasers  that  have  dif¬ 
ferent  angular  divergences  in  two  dimensions.  One 
could  also  imagine  using  the  spatially  dependent  light 
shift  of  atoms  traveling  through  the  ring  pattern  as  a 
Fresnel  lens  for  focusing  the  atoms.  One  future  topic 
for  study  is  the  long-distance  propagation  of  the  sub- 
Gaussian  transmitted  beam  and  the  super-Gaussian 
reflected  beam.  We  have  calculated  that  the  phase 
fronts  of  the  flat  profile  region  are  quite  uniform  and 
thus  expect  this  portion  of  the  beam  to  propagate  as 
a  plane  wave  over  a  long  distance. 

ThiB  research  was  supported  by  the  National  Sci¬ 
ence  Foundation,  the  U.S.  Office  of  Naval  Research, 
and  the  U.S.  Air  Force  Office  of  Scientific  Research. 
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Traisieflt  Laser  Cooling 
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We  observed  a  new  type  of  sub- Doppler  cooling  that  employs  neither  polarization  gradients  nor  mag.* 
netic  fields,  and  involves  neither  a  damping  force  nor  significant  diffusive  heating.  Instead,  light  shifts 
combined  with  optical  pumping  (OP)  to  levels  not  coupled  by  the  laser  field  gives  transient  cooling,  and 
steady  state  is  not  achieved.  The  time  scale  is  set  by  the  OP  rate.  We  observe  both  cooling  and  heating, 
but  often  with  the  opposite  detuning  from  that  of  steady  state.  A  semklassical  and  a  fully  quantum 
mechanical  calculation  of  transient  cooling  agree  very  well  with  one  another  and  with  our  data. 

PACS  nun  ben:  32.IO.Pj 


Laser  cooling  of  neutral  atoms  requires  velocity- 
dependent  optical  forces  that  arise  from  the  motion  of  an 
atom  in  a  nearly  resonant  light  field.  If  the  selection 
rules  require  that  excited  atoms  spontaneously  decay  to 
only  the  initial  state,  the  velocity-dependent  interaction 
arises  from  the  Doppler  shift  of  the  laser-driven  transi¬ 
tion.  In  this  case  one  readily  calculates  a  minimum 
achievable  temperature  called  the  Doppler  limit.  To 
where  r*l/y  is  the  excited-state  lifetime  III. 
By  contrast,  if  the  ground  state  has  multiple  sublevels  ac¬ 
cessible  from  the  excited  state,  much  lower  temperatures 
are  achievable  (2-4].  The  transitions  from  these  sublev¬ 
els  may  be  driven  by  different  polarizations,  and  thus 
sub-Doppler  temperatures  are  often  associated  with  po¬ 
larization  gradients.  The  cooling  force  derives  from  the 
failure  of  the  internal  state  of  moving  atoms  to  follow 
adiabatically  the  changing  optical  environment  of  a  spa¬ 
tially  inhomogeneous  optical  field,  such  as  that  produced 
by  light  beams  of  different  polarizations. 

For  slow  enough  atomic  velocities,  however,  this  non- 
adiabatic  part  of  the  force  is  small  compared  with  the 
velocity-independent  force  that  derives  from  the  spatially 
varying  light  shifts  of  the  atoms.  In  a  standing  wave,  this 
conservative  force  produces  channels  that  can  strongly 
influence  atomic  motion  or  even  confine  atoms  in  X/2  size 
regions  ($.6).  Furthermore,  spontaneous  decay  in  a 
standing  wave  can  provide  the  irreversible  process  needed 
to  produce  a  dissipative  force  that  cools  atoms  (7]. 

We  report  here  a  theoretical  and  experimental  study  of 
tub- Doppler  atom  cooling  achieved  by  transient  effects  of 
laser  excitation.  We  have  found  that  cooling  and  heating 
can  occur  when  the  atom  goes  through  only  a  few  opti¬ 
cal  pumping  (OP)  cycles,  but  often  with  opposite  de¬ 
tuning  from  that  of  steady-state  situations.  In  most  other 
types  of  sub-Doppler  laser  cooling,  a  damping  force 
F*-tr(pfW)  can  be  computed  from  the  steady-state 
solution  of  the  optical  Bloch  equations  for  the  density  ma¬ 
trix  p.  This  is  appropriate  because  the  laser  cooling  pro¬ 
cesses  generally  continue  for  times  long  compared  with 
OP  transients.  In  transient  laser  cooling,  the  time  scale  is 
set  by  the  OP  rate  to  a  ground  state  that  is  not  coupled  to 
an  excited  state  by  the  laser  field  (“dark  state")  to 
steady-state  solutions  are  inappropriate  and  a  damping 


force  cannot  be  defined  (8]. 

For  an  intuitive  understanding  of  how  transient  cooling 
works,  consider  the  conservative  motion  of  atoms  that 
enter  a  standing-wave  field  nearly  transverse  to  its  k  vec¬ 
tors,  but  with  a  small  velocity  component  r,  of  a  few  cm/s 
parallel  to  k  “XL  Such  atoms  experience  a  spatially 
varying  light  shift  potential  t/(*)  •»C/osin  Hks)  that  pro¬ 
duces  a  sinusoidal  force.  Wo  takes  the  sign  of  the  detun¬ 
ing  dwwutT.-a^—  and  is  proportional  to  the  intensity  at 
low  intensity.) 

The  atoms  experience  this  force  until  they  are  optically 
pumped  to  a  dark  state.  Since  the  OP  rate  y,(z)»y, 
x  sin  JUz)  vanishes  at  the  nodes,  OP  is  slower  for  atoms 
that  enter  the  standing  wave  near  a  node,  and  thus  the 
average  kinetic  energy  (KE)  change  is  greater  for  such 
atoms.  For  S>  0,  more  KE  is  lost  by  atoms  moving  away 
from  a  node  than  gained  by  atoms  moving  away  from  an 
antinode,  so  on  average  atoms  lose  KE  (Fig.  I  (a)]  and 
conversely  for  6  <  0  (Fig.  1(b)].  This  contrasts  with  oth¬ 
er  types  of  cooling  where  the  optical  force  may  either 
heat  or  cool  the  atoms  depending  on  both  the  angular 
momentum  scheme  (Ff-fV)  end  the  detuning:  Here 
S  >  0  always  gives  cooling. 

We  have  performed  classical  trajectory  calculations 
based  on  the  above  model  by  integrating  the  equation 
mi' — k(/o$in(2kz).  For  initial  conditions  we  choose  i<o 
and  zo  distributed  uniformly  in  the  intervals  ("30,30] 
cm/s  and  (0, X/2],  respectively.  In  the  weak  excitation  ap¬ 
proximation,  the  upper  state  can  be  neglected,  and  the 
ground-state  potential  amplitude  Uom2HsCS/L ,  where 
sml/!m  is  the  laser  saturation  parameter,  C  is  the 
strength  of  a  particular  transition,  and  £.•!+ (25/y)1. 
The  saturation  intensity  for  the  strongest  component  of 
this  transition  (C"l)  b  /^•kf/X’r  a  1.6  mW/cm*  in 
Rb,  where  t  -  l/y»  27  ns. 

From  the  solution  sit)  to  these  equations  of  motion  we 
find  the  OP  rate  y,(f)«f#sin*ttz(r)],  where  y,m2sC 
xByliL+sC)  and  £  is  the  branching  ratio  to  the  dark 
Mates.  Then  the  probability  of  OP  between  l  and  r+rfi  is 
Pit)dt - y,(r)expt -/4y, ((')</« ']di.  The  primitive  of 
Pit)  is  then  inverted  to  find  a  distribution  of  pumping 
times.  A  Monte  Carlq  approach  uses  these  pumping 
times  to  terminate  the  classical  motion,  and  the  resulting 
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where  OP  ocean,  we  use  dt  mdz/v,  end  set 


FIG.  I.  (a)  The  tight  shift  potential  and  optical  pumping 
scheme  for  i>  0.  Atoms  entering  near  a  node  travel  further 
before  OP,  and  therefore  loae  more  KE  than  atoms  entering 
near  an  antinode,  which  travel  a  shorter  distance  before  OP. 
(b)  Similar  to  (a),  but  for  S  <0  atoms  gain  more  KE  than  they 
lose,  (c)  Velocity  distributions  found  from  the  aemiclassica! 
calculation  (dots)  for  the  transition  shown  in  the  inset.  Here 
i«0.66  and  8 •  + 12  MHz  (as  in  Fig.  3).  “Snapshots"  of  the 
velocity  distribution  are  made  at  times  3, 6,  and  9  ps.  The  solid 
lines  show  the  results  of  quantum  mechanical  calculations  for 
the  same  conditions,  (d)  Similar  to  (c)  but  for  6*  - 12  MHz. 
Inset:  One  of  the  simplest  level  schemes  for  transient  cooling, 
Ft  m  I  -*  Ft  m0  with  a*  excitation. 

velocity  distribution  it  calculated. 

We  have  applied  this  approach  to  thn  simplest  atomic 
scheme  of  the  classical  scheme  above,  an  Fm\  -*  0  tran¬ 
sition  with  a*  excitation,  shown  in  Fig.  I  (inset).  Re¬ 
sults  of  these  calculations  are  shown  in  Figs.  1(c)  and 
1(d)  (dots).  The  sharp  features  in  Fig.  1(d)  for  5  <  0 
occur  at  stationary  points  where  £**  **(/<>.  It  is  clear 
that  this  simple  model  of  transient  User  interactions  gives 
cooling  for  blue  detuning  and  heating  for  red  detuning  for 
this  £»  I  -*  0  transition. 

Figures  I  (c)  and  I  (d)  also  show  remarkable  agreement 
between  the  class*  il  results  and  results  of  a  fully  quan¬ 
tum  mechanical  calculation  including  the  excited  state 
explicitly  (solid  lines).  The  sharp  cooling  and  heating 
features  in  the  classical  results  are  smoothed  in  the  quan¬ 
tum  results.  We  attribute  this  to  neglect  of  the  uncer¬ 
tainty  in  position  and  of  recoil  in  the  classical  calculation. 

We  can  estimate  the  velocity  capture  range  of  transient 
blue  cooling  for  atoms  that  move  freely  without  being 
channeled  between  the  planes  of  the  standing  waves.  The 
probability  of  OP  in  time  »i  is  F(t\)mfi'Y,U)dt  for 
P(fi)«l.  To  find  the  mean  distance  xi  from  a  node 


J*'l  yfiz)ht\dz  -  • 

For  E  >  21/©,  the  variation  of  r,  with  z  has  little  effect  on 
z  i  so  it  comes  out  of  the  integrals,  and  we  find 
sin  :Uri)  “0.83.  Thus  atoms  undergo  OP  from  a  posi¬ 
tion  where  their  average  light  shift  (potential  energy)  is 
—  f  Vo,  and  since  their  initial  spatial  distribution  was 
uniform  so  their  average  initial  potential  energy  was 
VJ2,  the  average  total  energy  change  is  A£ — t/o/3. 
This  is  a  mechanical  energy  loss  for  l/o  >  0  (S  >  0)  and  a 
gain  for  t/o<0(£<0).  This  process  does  not  make  a 
very  large  change  in  the  velocity  distribution  for  atoms 
with  E  more  than  a  few  times  (/©  because  Ar,/vt  “A £/ 

2 E - CV6£.  This  estimate  of  the  capture  velocity  is 

consistent  with  the  detailed  computational  results  of  Fig. 
I  and  with  the  data. 

In  this  transient  cooling  process,  the  final  velocity  dis¬ 
tribution  does  not  result  from  competition  between  a 
steady-state  damping  force  and  diffusive  heating.  Instead 
we  find  the  changes  in  KE  are  bounded  by  t/o.  The 
widths  of  the  cooling  peaks  (and  heating  dips)  in  the 
measured  velocity  distributions  decrease  with  intensity  so 
that  widths  substantially  below  the  Doppler  limit  are  at¬ 
tainable,  as  shown  by  our  measurements. 

The  experiments  use  a  thermal  beam  of  natural  Rb 
produced  by  an  oven  at  7*—  150°C  with  horizontal  slit 
aperture  0.06  mm  high  by  2  mm  wide,  and  a  vertical 
beam  defining  slit  2  mm  high  by  0.06  mm  wide  about  35 
cm  away  (9].  The  atoms  emerge  from  the  vertical  slit  in 
a  fan-shaped  beam  and  then  interact  with  a  pair  of  coun- 
terpropagating  laser  beams  transverse  to  the  atomic  beam 
axis.  The  nearly  flat  atomic  beam  profile  is  measured 
with  a  scanning  hot  platinum-tungsten  wire,  25  pm  in  di¬ 
ameter,  1.3  m  away  from  the  region  of  interaction  with 
the  laser  beam.  Three  square  Helmholtz  coil  pairs  cancel 
the  Earth’s  field. 

A  35  mW  Sharp  LT02S  diode  laser  is  locked  at  the 
5 S—  SP  transition  of  Rb  near  X«780  nm  by  saturated 
absorption  in  a  vapor  cell  at  room  temperature.  The  spa¬ 
tial  intensity  profile  of  the  23  mm  long  *5  mm  high  laser 
beam  is  flattened  to  a  few  percent  by  a  tilted  ctalon  llOl. 
For  thermal  velocity  atoms  (F— 350  m/s)  the  average  in¬ 
teraction  time  is  63  gas.  The  tum-on  and  turn-off  edges  of 
the  light  field  are  sharpened  by  an  aperture  and 
diffraction  limited  to  about  0.3  mm,  corresponding  to  an 
eatranoe  and  exit  time  of  <  I  its  for  atoms  at  thermal  ve¬ 
locity.  For  most  experimental  conditions,  the  OP  time  is 
longer  than  I  pa 

In  order  to  explore  the  phenomena  sketched  above,  we 
have  measured  the  transverse  velocity  distribution  of  our 
atomic  beam  after  it  passes  through  the  circularly  polar¬ 
ized  aunding  wave  in  zero  magnetic  field.  For  the 
Ft  m  I  —  Ft  *2  transition  of  *’Rb  (dashed  line  in  Fig.  2), 
OP  to  the  dark  sute  Ft  “2  is  allowed.  Figure  3(a)  shows 
clear  evidence  of  cooling  for  light  detuned  blue  by 
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£-  +  12  MHz  (»2y)  from  this  transition,  and  Fig.  3(b) 
•hows  heating  for  ff-  — 12  MHz.  This  demonstration  of 
blue  cooling  (sub-Doppler)  and  red  heating  appears  in  all 
our  data  for  all  the  allowed  transitions  of  Fig.  2. 

Figure  3(a)  shows  significant  blue  cooling  even  with  s 
lowered  to  —0.066  by  a  neutral  density  (ND— 2)  filter. 
At  this  intensity,  the  average  OP  time  f,"  l/y>,  found  by 
averaging  yp  over  a  wavelength  and  over  the  various  tran¬ 
sitions,  is  64  /is.  (Typically,  yrlyf  ranges  from  2  to  4.) 
This  is  fortuitously  close  to  the  average  63  /is  interaction 
time  for  atoms  in  our  thermal  beam.  When  the  interac¬ 
tion  region  is  shortened  from  the  23  mm  of  Fig.  3(a)  to 
only  1.4  mm,  corresponding  to  an  average  interaction 
time  of  only  4  /is,  there  is  still  dear  evidence  of  cooling 
on  the  blue  side  for  s*6.6. 

Excitation  of  ,5Rb  on  the  Ff-2— >F,-1  transition 
(solid  line  of  Fig.  2)  presents  a  dose  analog  to  the  Ft 
-  I  -*  Ft  «0  model  discussed  above,  because  the  F,  - 1 
state  cannot  decay  to  Ff-  3.  However,  OP  can  occur  to 
the  F,- 2,  Mf-  +  I  and  +2  magnetic  sublevels,  which 
•re  “dark"  in  our  a*  light.  The  low  intensity  data  near 
the  top  of  Fig.  4(a)  clearly  show  cooling  signals  with  <r4 
light  detuned  +6  MHz  blue  of  this  F, -2—  F,  -  1  tran¬ 
sition.  At  the  lowest  intensity  of  Fig.  4(a)  (top  trace), 
l/o lh  is  only  1 1  kHz  so  the  expected  change  in  the  veloci¬ 
ty  of  an  atom  is  less  than  the  recoil  velocity  hk/M,  the 
transient  cooling  effects  are  washed  out  by  the  momen¬ 
tum  changes  from  the  random  direction  of  spontaneous 
emission. 

At  higher  intensity,  a  second  phenomenon  appears. 
Light  tuned  6  MHz  blue  of  this  Ftm2—  F,-l  transi¬ 
tion  of  MRb  is  also  23  MHz  red  of  the  F#"2—  F,m 2 
transition  (dotted  line  of  Fig.  2).  Although  atoms  opti¬ 
cally  pumped  to  the  (2,1)  sublevel  cannot  be 

excited  \oF,m\  with  a*  light,  there  is  an  allowed  transi¬ 
tion  to  the  (F#,Mr)—(2,2)  sublevel.  Since  this  state  can 
decay  to  the  F,m  3  dark  state  (wiggly  line  in  Fig,  2), 
transient  effects  can  occur.  Using  rf  calculated  as  be¬ 
fore,  we  find  the  rate  for  this  further  off-reaonant  transi¬ 
tion  to  be  significant  only  at  high  enough  intensity.  This 
slower  process  involves  light  tuned  red  of  resonance,  so  it 
heats  instead  of  cools  the  atoms  as  shown  in  the  high  in¬ 
tensity  traces  of  Fig.  4(a).  The  data  of  Fig.  4(a)  clearly 
show  a  change  from  cooling  to  heating  as  a  increases 


FIG.  3.  (a)  Measured  atomic  beam  profiles  for  circularly  po¬ 
larized  light  tuned  12  MHz  blue  of  the  Ff-l  —  F,-2  transi¬ 
tion  of  rRb.  The  intensities  are  0.1,  0.3,  1.0,  3.0,  and  10 
mW/cm2  (bottom  to  top),  and  the  saturation  intensities  for  the 
three  magnetic  components  are  19.2, 6.4,  and  3.2  mW/cm2  (3.2 
for  Mr— 1  —  2).  The  vertical  scale's  arbitrary  units  are  — 1% 
of  the  atomic  beam  intensity.  For  the  mean  longitudinal  veloci¬ 
ty  of  3S0  m/s,  a  37  /im  displacement  represents  I  cm/s,  so  the 
width  of  the  third  peak  (saturation  *1  for  M-l  —  2)  is  6 
cm/s,  corresponding  to  about  i  of  the  Doppler  limit,  (b)  Simi¬ 
lar  reiults  for  light  tuned  12  MHz  ted  of  the  same  transition. 


from  0.3  to  1.0  for  our  63  /is  interaction  time,  corre¬ 
sponding  to  tp  going  from  80  to  24  /is. 

We  can  also  observe  this  phenomenon  by  varying  the 
interaction  time  at  fixed  s.  Figure  4(b)  shows  transient 
oooling  on  this  Ftm2-*Ftm\  transition  at  short  times. 


FIG.  4.  (a)  Measured  atomic  beam  profiles  for  circularly  po¬ 
larized  light  tuned  6  MHz  blue  of  the  F,-2—F, -I  transition 
of  “Rb.  The  iateMKias  are  10.  3.0,  1.0.  0.3,  0.16,  0.03,  and 
0.01  mW/cm1  (bottom  to  top),  and  tbe  saturation  intensities  for 
the  three  components  are  16.  S.I.  tad  2.7  mW/cm 2  (16  for 
Mr-0—  I).  For  the  mean  longitudinal  velocity  of  350  m/s,  a 
37  pm  displacement  represents  I  cm/s.  so  the  width  of  the 
fourth  peak  ie  4  cm/s,  corresponding  to  about  t  of  the  Doppler 
limit  (b)  Similar  data  for  interaction  timet  of  32, 16, 8,  and  4 
pi  (bottom  to  top).  Cooling  becomes  beating  at  higher  intensi¬ 
ty  or  longer  interaction  time. 
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FIG.  3.  Velocity  distributions  found  from  the  fully  quantum 
calculation  with  F, -2,  F# -3,  F,-l,  and  F,  -2  internal  states 
included,  using  appropriate  excitation  and  decay  m.  fix  ele¬ 
ments.  The  saturation  parameter  is  taken  to  be  0.4  and  the  de- 
tun'ng  is  6  MHz  above  the  F, “I  state,  as  for  the  experimental 
data  in  Fig.  4(b).  The  successive  curves  correspond  to  interac¬ 
tion  umcs  of  4, 8, 16,  and  32  /is  for  the  most  probable  longitudi¬ 
nal  velocity  in  the  beam. 


but  heating  as  the  interaction  time  increases  from  16  to 
32  /is  at  s  •0.35.  In  this  case  atoms  in  2,1) 

are  given  enough  time  for  excitation  to  F,  -2.  There  is 
also  a  small  but  persistent  and  repeatable  upward  peak 
from  the  bottom  of  the  heating  dip  in  the  high  intensity 
data  of  Fig.  4(a)  and  in  the  long  time  data  of  Fig.  4(b). 
It  may  come  from  the  (Ff,A#f)-(2,2)  atoms  that  are 
not  heated  by  excitation  to  Ft  -2.  All  the  features  of 
Fig.  4  have  also  been  observed  on  the  F,  ■  1  —  F,  “0 
transition  of  ,7Rb. 

As  a  further  illustration  of  this  transient  cooling  model, 
we  have  added  a  small  magnetic  field  transverse  to  the 
axis  defined  by  the  circularly  polarized  standing  wave 
(BXk).  When  the  Larmor  frequency  is  large  enough  to 
precess  atoms  in  the  {FtMr )  -(2, 1 )  or  (2,2)  dark  states 
back  to  the  coupled  (2,0)  state  faster  than  they  can  be 
optically  pumped  to  the  Ff  »3  dark  state  (via  F,  —2,  wig- 
gly  line  in  Fig.  2),  the  cooling  persists  to  the  higher  inten¬ 
sities  shown  in  Fig.  3(a). 

Multistate  processes  are  difficult  to  model  accurately 
with  classical  trajectory  calculations,  but  we  have 
developed  a  fully  quantum  mechanical  model  (I  I).  The 
change  from  cooling  to  heating  shown  in  Fig.  4  provides 
an  opportunity  for  sensitive  tests  of  such  quantum  me¬ 
chanical  calculations.  The  basis  set  consists  of  product 
states  of  internal  and  external  atomic  coordinates  (12]. 
The  set  of  internal  states  includes  typically  two  ground 
and  two  excited  levels  so  the  calculation  is  not  restricted 
bo  low  excitation  rates.  The  external,  center-of-mass 
motion  states  are  free  particle  wave  functions,  allowing 
inclusion  of  the  recoil  effect  in  every  optical  transition. 
The  density  matrix  calculated  in  this  basis  includes  opti¬ 
cal  coherence  (off-diagonal  elements)  to  sufficiently  high 
order  so  that  the  calculated  velocity  distribution  con¬ 


verges.  Equations  for  the  temporal  evolution  of  the  densi¬ 
ty  matrix  are  obtained  from  the  optical  Bloch  equations. 

Figure  3  shows  the  results  from  these  quantum  calcula¬ 
tions  for  parameters  corresponding  to  Fig.  4(b).  The  ex¬ 
perimental  conditions,  including  the  longitudinal  velocity 
distribution  and  instrumental  resolution,  have  been  care¬ 
fully  modeled.  The  transverse  velocity  distribution  flips 
from  cooling  to  heating  as  in  the  experiments,  but  there  is 
an  unexplained  discrepancy  of  about  a  factor  of  2  in  the 
*.-»h  of  the  final  heating  dip. 

There  are  various  repumping  methods  that  might  allow 
recycling  of  atoms  for  further  cooling  by  this  technique 
(7).  Furthermore,  it  can  readily  be  extended  to  2D  or 
maybe  even  3D.  Also,  this  technique  might  be  applied  to 
atoms  with  a  first  transition  in  the  deep  ultraviolet  and  a 
metastable  state  where  traditional  cooling  techniques  are 
difficult  to  implement.  For  example,  atomic  hydrogen  in 
its  2 S  state  could  be  cooled  this  way  on  the  easily  pro¬ 
duced  Balmer-c  line  at  656  nm.  Similar  possibilities  exist 
in  singlet  He,  rare  gases,  and  alkaline  earth  D  states. 

In  summary,  our  experiments  with  a  Rb  atomic  beam 
and  a  transverse  standing-wave  laser  field  have  demon¬ 
strated  a  new  mechanism  for  sub-Doppler  laser  cooling. 
This  mechanism  requires  no  polarization  gradients  or 
magnetic  fields.  It  operates  on  a  short  time  scale,  deter¬ 
mined  by  the  OP  times.  We  have  developed  a  theoretical 
model  with  ciassical  trajectories  and  have  also  applied  a 
fully  quantum  mechanical  treatment.  These  models  give 
results  close  to  each  other  and  to  the  measurements. 

This  work  was  supported  by  NSF,  ONR,  AFOSR,  and 
CAPES  (Brazil).  The  quantum  calculations  were  per¬ 
formed  at  the  Cornell  National  Supercomputer  Facility, 
funded  by  NSF  and  IBM. 
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Quantum  calculations  for  ID  laser  cooling: 
temporal  evolution 

T.  Bergeman 

Department  of  Physics,  State  University  of  New  York 
Stony  Brook,  NY  11194-3800 


Temporal  integration  of  density  matrix  equations  with  quantized  translational 
motion  gives  new  results  for  laser  cooling  of  atoms.  For  polarization  gradient  cool¬ 
ing,  the  average  kinetic  energy  in  the  cold  atom  peak  after  1,000  radiative  lifetimes 
decreases  with  increasing  Fg,  and  is  <1  recoil  energy  at  low  laser  intensity.  With 
magnetically  induced  laser  cooling  (M1LC)  and  low  optical  pumping  rates,  atoms 
cool  slowly  within  the  potential  wells,  giving  minima  in  the  velocity  distribution  at 
v=0.  MILC  is  seen  as  a  cyclic  process;  optical  pumping  cools  and  B  field  mixing 
slightly  heats. 


In  view  of  many  recent  advances  in  techniques  for  cooling  atoms  by  laser  light, 
there  is  interest  in  theoretical  methods  that  can  address  diverse  experimental  situa¬ 
tions.  Most  of  the  theoretical  results  presented  to  date  have  been  for  the  steady-state 
limit  of  the  cooling  process.  For  example,  semidassical  theories(1.2]  for  atom  cool¬ 
ing  employ  a  force  versus  velocity  function  and  a  diffusion  parameter  based  on  the 
steady-state  solution  of  the  density  matrix  equations.  The  first  examples  of  density 
matrix  calculations  with  a  basis  of  quantized  translational  states  [3.4]  also  presented 


steady-state  solutions,  obtained  by  inverting  the  evolution  matrix.  Exceptions  to 
the  steady-state  limit  are  offered  by  quantum  Monte  Carlo  methods,  which  have 
successfully  reproduced  laser  heterodyne  results  [5].  However,  these  methods  have 
not  been  extensively  used  for  multi-level  atomic  transitions  because  of  the  large 
number  of  replications  needed. 

In  many  ID  atom  cooling  experiments  where  an  atomic  beam  traverses  coun- 
terpropagating  laser  beams,  steady-state  solutions  are  not  suitable  for  comparisons 
with  experiment.  For  atoms  that  have  been  precooled  in  a  3D  magneto-optic  trap, 
the  interaction  times  in  a  ID  standing  wave  may  indeed  be  adequate  to  achieve 
steady  state  [6].  However,  in  experiments  with  a  thermal  atomic  beam,  the  inter¬ 
action  time  is  typically  no  more  than  1,000  to  2,000  r,  where  r  is  the  excited  state 
radiative  lifetime.  The  steady  state  limit  may  not  be  attained,  particularly  when 
the  initial  velocity  distribution  is  appreciably  wider  than  the  so-called  capture  ve¬ 
locity,  vc.  The  range  of  atom  velocities  that  are  effectively  slowed  depends  on  the 
interaction  time.  Clearly  a  theory  is  needed  for  finite  interaction  times. 

It  will  be  demonstrated  in  this  report  that  density  matrix  methods  over  a  basis 

t 

of  free  particle  states  or  of  eigenfunctions  of  the  light  shift  potential  are  capable 
of  significant  extension  beyond  their  use  as  previously  reported  [3. 4).  By  efficient 
programing  and  the  use  of  supercomputer  technology,  results  have  been  obtained 
for  the  temporal  evolution  and  for  higher  angular  momentum.  The  computational 
techniques  will  be  illustrated  by  both  lin  1  lin  polarization  gradient  cooling  and 
magnetically  induced  laser  cooling  ( MILC).  Results  will  be  shown  for  cases  in  winch 
the  incremental  velocity  change  (the  recoil  velocity)  is  not  small  compared  with 
the  width  of  the  velocity  distribution,  thus  where  semiclassical  Fokker-Flanck  ap¬ 
proaches  would  not  be  valid.  Quantum  methods  will  also  be  applied  to  atoms 
whose  kinetic  energy*  is  less  than  the  amplitude  of  the  periodic  light  shift  potential. 


In  this  regime,  the  use  of  a  constant  velocity  parameter  in  the  computation  of  the 
semiclassica!  force  function  is  inappropriate. 

Calculations  with  a  free  particle  basis  here  include  the  excited  state  explicitly 
(3].  The  basis  states,  |F,  mj?,  n  >=  |F,mf  >  txp(inks/h ),  where  2n/k  is  the  laser 
wavelength,  are  products  of  internal  atomic  sublevels  and  free-particle  momentum 
eigenfunctions.  In  the  evolution  equation  for  the  density  matrix,  a  —  —(i/h)[H,  <r]+ 
<rrt|,  the  Hamiltonian,  H,  includes  the  laser-atom  interactions  (Hl)  and  possibly 
magnetic  field  effects,  while  <r,eJ  expresses  the  effects  of  spontaneous  decay  and 
repopulation,  including  an  average  over  the  spatial  distribution  of  emitted  light. 
Elements  of  Hl  include  photon  recoil  explicitly.  Elements  of  a  are  placed  in  a 
column  vector,  x,  and  the  coupled  equations  x  =  W’x  are  integrated  numerically. 
These  computer  programs  have  been  modified  to  deal  with  transitions  with  AF 
=  -1  and  0,  more  than  one  ground  and  excited  state  hvperfine  level,  and  spatially 
nonuniforxn  laser  intensity.  In  this  report,  however,  only  F  — »  F  -f  1  transitions  are 
considered  and  the  laser  intensity  is  constant  over  the  interaction  region. 

A  full  density  matrix  for  up  to  ±50  recoil  momenta,  even  for  a  mesh  spacing 
of  one  recoil  and  a  7  =  1/2  -*  3/2  transition,  would  have  more  than  3  x  10* 
elements  and  the  evolution  matrix  would  then  have  10s*  elements.  However,  off’ 
diagonal  elements  with  An  >  n„er  are  found  not  to  affect  the  computed  velocity 
distribution.  varies  with  laser  intensity  and  detuning  but  is  typically  6  to 

12.  Effectively,  a  becomes  a  band  matrix  of  order  n.  In  the  present  calculations, 
there  were  up  to  100,000  density  matrix  elements.  The  evolution  matrix.  W.  can 
be  simplified  by  ueglccting  elements  less  than  10'1J  in  magnitude,  for  example. 
No  more  than  20  to  50  elements  are  needed  in  each  row  for  the  free  particle  basis 
approach.  By  such  means,  it  is  possible  to  compute  the  teui|>oral  evolution  to  1.000 
radiative  lifetimes  with  a  basis  of  65  momentum  values  for  each  F.  my  sublevel  in 
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an  F  =  3  — *  4  transition  with  about  2  hours  of  equivalent  Cray  XMP  CPU  time. 

When  the  atoms  are  weakly  excited,  excited  state  elements  may  be  eliminated 
from  the  density  matrix  equations  to  obtain  a  periodic  light  shift  potential  for  the 
ground  state  [4].  This  approach  gives  a  useful  physical  picture,  and  for  certain 
applications,  uses  less  computer  time.  The  energy  spectrum  exhibits  a  band  struc¬ 
ture  resembling  that  obtained  for  electron  waves  in  a  ID  crystal  lattice.  Optical 
pumping,  magnetic  field  mixing,  and  spontaneous  emission  rates  are  computed  for 
density  matrix  elements  over  the  basis  set  of  Bloch  states. 

Results  obtained  for  the  lin  1  lin  laser  configuration,  in  which  two  counterprop- 

agating  laser  beams  have  orthogonal  linear  polarization,  show  that  the  steady-state 

limit  may  not  reliably  represent  cooling  processes  occuring  over  shorter  interaction 

times.  Fig.  1  shows  typical  time  evolution  plots  obtained  with  the  free  particle 

basis.  Here,  e  =  E^/hT  —  6.4  x  10“*.  as  for  Rb.  where  Er  =  h2k2/2M  is  the  recoil 

energy  and  F  is  the  radiative  decay  rate.  The  laser  field  is  characterized  by  the 

light  shift  well  depth.  Uq  —  - fShi/L .  For  lin  1  lin  cooling  on  an  F  — *  F  +  1 

transition.  /  =  [(2F+  1)(F+  1)  -  1]/[(2F+  1)(F  +  1)},  5  =  2Q*/T2  where  Q  is  the 

* 

single  beam  Rabi  frequency,  b  =  wj#M,  -  is  the  detuning,  and  1=1  +  462/F2. 
(5  =  1  for  laser  intensity  I  =  hc/A3r).  The  calculations  shown  in  Fig  la.  for 
F  =  1/2  —  3/2.  are  carried  to  100  pset\.  or  3770  r.  where  r  =  1/r.  Ultimately, 
the  velocity  distribution.  Pit),  does  evolve  into  the  typical  form  for  the  asymptotic 
hunt, namely  a  narrow  Gaussian  on  top  of  a  broad  Gaussian  (3j .  However,  up  to  40 
ps.  the  narrow  Gaussian  is  flanked  by  two  dips  at  tlie  capture  velocity,  er  (about 
16  v r  here).  The  time  required  to  attain  the  steady  state  distribution  depends  on 
the  width  of  the  initial  distribution  (here  assumed  flat  up  to  »  CO)  as  well  as  on 
the  intensity  and  detuning. 

With  lower  laser  intensity  or  wider  initial  distribution  over  t\  the  dip*  at  r  at 


±vc  persist  for  longer  times.  The  generic  behavior  is  shown  in  Fig.  lb,  for  F  = 
3  — »  4.  To  extract  an  average  kinetic  energy  for  atoms  in  the  central  peak,  two 
procedures  have  been  used.  The  first  assumes  that  the  distribution  for  |u|  <  vc  is 
Gaussian.  The  velocity  vx/t  at  1/e  of  the  maximum  height  difference  then  yields  a 
Gaussian  average  energy,  KEq  =  Mvl/J2.  In  the  second  method,  one  computes  the 
average  kinetic  energy,  KEa,  of  atoms  in  the  range  |u|  <  vc-  When  the  distribution 
in  the  central  peak  is  truly  Gaussian,  KEa  —  KEg,  but  typically,  KEa  ~  O.lKEc- 

The  variation  of  KEq  with  laser  parameters  after  30  /is  (1,130  r  for  Rb)  is 
shown  in  Figs.  2a  and  2b.  For  the  calculations  shown  in  Fig.  2,  e  =  6.4  x  10-4, 
6  =  —  5F,  and  again  the  free  particle  basis  was  used.  To  emphasize  the  role  of  the 
well  depth,  in  Fig.  2a  I\Eg/U0  is  plotted  vs.  Uq/Er.  For  F3  =  3,  the  average 
kinetic  energy  in  the  cold  atom  peak  is  significantly  less  than  for  Fg  =  1/2.  Since 
this  trend  seems  to  persist  to  longer  interaction  times,  this  may  in  part  explain  the 
experimental  observation  that  the  fraction  of  atoms  in  the  lowest  quantum  state 
in  experiments  on  ^Rb,  with  Fg  =  3,  was  larger  than  computed  for  Fg  —  1/2  [6]. 
One  possible  explanation  for  the  lower  kinetic  energies  with  high  Fg  is  the  greater 

t 

importance  of  Amp  =  ±2  coherences  induced  by  the  laser  field,  which  produce 
progressively  larger  off-diagonal  elements  in  the  light-shift  potential  matrix.  For 
high  F,  energy  bands  within  the  potential  well  are  broadened,  making  the  motion 
for  I\E  <  Uq  more  nearly  like  that  of  free  particles. 

For  finite  interaction  times  and  Uq/Er  <  100,  the  average  kinetic  energy  for 
atoms  in  the  central  peak  is  significantly  below  that  indicated  iu  steady-state  results, 
which  gave  a  minimum  of  KE  5 Er  (4).  Here,  there  appears  to  be  no  minimum, 
and  at  low  intensity  (Fig.  2b)  KEa  <  Er.  Because  the  amplitude  of  the  cooling 
peak  is  only  5  •  10%  of  the  initial  amplitude  at  v=0  for  tins  regime  of  low  Uq/Er. 
such  extremely  cold  atom  distributions  may  be  difficult  to  observe. 
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For  MILC  [7,  8],  the  light  shift  potential  matrix  is  diagonal  and  sinusoidal, 
so  the  wave  equation  for  each  mf  is  simply  a  Mathieu  equation.  The  periodic 
potential  basis  approach  is  appropriate  for  S/L  <<  1.  To  model  experiments  in 
which  atoms  move  freely  outside  the  interaction  region,  an  initially  fiat  distribution 
over  free  particle  states  is  projected  onto  the  periodic  state  basis.  After  computing 
the  evolution  to  times  the  distribution  over  the  periodic  state  basis  is  projected 
back  onto  the  free  particle  basis.  Tr(<r)  decreases  by  a  few  percent  over  typical 
interaction  times  because  the  total  optical  pumping  rate  out  of  any  quantum  basis 
state  is  greater  than  the  repopulation  terms  from  this  level  summed  over  any  finite 
set  of  basis  states.  Results  from  the  periodic  potential  basis  states  and  free  particle 
basis  states  are  very  similar  for  weak  excitation. 

The  temporal  behavior  computed  for  MILC  exhibits  a  new  feature  associated 
with  the  localization  of  atoms  in  the  standing  wave.  In  MILC,  the  eounterpropa- 
gating  laser  beams  have  identical  circular  polarization  and  the  light  shift  potential 
maxima  and  minima  for  different  mF  sublevels  are  spatially  in  phase.  A  transverse 
magnetic  field  mixes  different  m?  sublevels.  Zeeman  mixing  occurs  most  effectively 
at  the  nodes  of  the  laser  standing  wave,  where  the  potentials  are  degenerate.  (Mag¬ 
netic  field  mixing  serves  the  same  function  as  the  polarization  gradient  in  lin  1  lin 
cooling.)  With  red  detuning,  when  the  final  kinetic  energy  is  less  than  the  depth 
of  the  light  shift  potential  wells  in  MILC,  atoms  in  the  ground  state  cannot  tra¬ 
verse  the  potential  peaks  at  the  nodes,  the  magnetic  field  mixing  is  suppressed,  and 
the  cooling  rate  slows.  Although  the  asymptotic  velocity  distribution  may  still  be 
nearly  Gaussian,  over  very  extended  times  there  is  a  deficiency  of  the  coldest  atoms. 
The  velocity  distribution  exhibits  a  flattened  top.  or  even  a  minimum  at  v  =  0.  Fig. 
3  shows  this  effect  in  the  velocity  distribution  for  F3  »  1/2  (in  Fig.  3a- c)  and  for 
Fi  s  3  (in  Fig.  3d).  In  MILC,  the  effective  optical  pumping  rate.  Fr  -  2 SV/L. 
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and  Uq  —  -2Sh6/L.  In  Fig.  3a  and  3c,  hTv  =  U0.  When  is  small  compared 
with  U0,  as  in  Fig.  3b  and  3d  for  which  hTp  =  Uo/20,  the  minimum  in  P(v)  at 
v  =  0  is  more  persistent.  Fig  4a  and  4b  show  distributions  over  the  states  in  the 
periodic  potential  corresponding  to  the  P(v)  distributions  shown  in  Figs.  3a  and 
3b.  The  energies  of  these  states  are  shown  at  the  top  of  Fig.  4.  Below  the  potential 
maxima,  U0/h  =  1  MHz  here,  there  are  quasi-discrete  states,  while  for  E  >  Uq,  the 
gaps  between  the  bands  are  narrow.  Although  cooling  increases  the  population  of 
the  states  lying  below  Uq,  the  population  of  the  lowest  states  increases  very  slowly, 
particularly  when  hTp  <  U0. 

For  stationary  states  in  the  light  shift  potential,  the  explanations  for  sub- 
Doppler  laser  cooling  derived  from  semiclassical  models  [2]  do  not  apply.  A  damping 
force  does  not  explicitly  enter  the  present  calculations.  In  the  quantum  represen¬ 
tation,  cooling  occurs  because  of  preferential  optical  pumping  out  of  higher-lying 
states  into  lower-lying  states.  The  way  this  works  in  MILC  is  shown  in  very  schemat¬ 
ically  in  Fig.  5  for  states  in  the  light  shift  potential.  Energy  bands  are  shown  in 

columns  for  m?  =  -1/2  (left  of  each  transition  pair)  and  for  mp  =  1/2  (right  of 

» 

each  pair).  The  magnetic  field  and  the  laser  interactions  redistribute  the  popu¬ 
lations  over  this  basis  set.  The  widths  of  the  arrows  indicate  the  relative  optical 
pumping  rates  or  B-field  mixing  coefficients  from  various  initial  states.  Not  shown 
is  the  much  slower  process  of  diffusive  heating  in  Amp  =  0  excitation  and  decay.  In 
the  optical  pumping  process,  on  average  the  emitted  photon  takes  away  more  energy 
than  is  added  by  the  absorbed  laser  photon.  This  is  due  in  part  to  the  difference  in 
the  nif  ss  ±1/2  potential  minima.  In  addition,  optical  pumping  favors  final  states 
localized  at  the  antinodes,  namely  the  lower-lying  states  for  red  detuning.  On  the 
other  hand,  the  magnetic  field  does  not  vary  spatially  over  the  interaction  region,  so 
the  redistribution  by  the  magnetic  field  is  not  so  sharply  peaked  at  An  =  0.  These 
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processes  constitute  a  cooling  cycle.  The  reduction  of  energy  in  the  Amp  =  1  opti¬ 
cal  pumping  part  of  the  cycle  is  greater  than  the  increase  of  energy  that  occurs  in 
AmF  =  -1  magnetic  field  mixing  and  in  Amp  =  0  diffusive  heating. 

Comparisons  of  the  results  present  above  with  experiment  are  obviously  of  great 
interest.  FYom  preliminary  results  [9],  it  appears  that  the  lin  -L  lin  results  are  in 
good  agreement  with  experiment  on  Rb  and  also  that  results  for  MILC  agree  well 
with  experiments  on  He  2  3S  [10].  However,  experimental  cooling  peaks  observed  for 
certain  MILC  experiments  with  Rb  [9]  are  narrower,  without  the  dip  at  v=0  found 
above.  One  possible  explanation  is  that  the  optical  wavefronts  in  the  apparatus  used 
for  Rb  are  not  ideal,  and  trapping  in  the  light  shift  wells  occurs  less  than  in  the 
calculations.  These  comparisons  and  comparisons  with  semiclassical  calculations 
will  be  discussed  in  future  publications. 

This  work  was  supported  by  NSF,  ONR,  and  AFOSR.  I  am  indebted  to  H. 
Metcalf  J.  Dalibard,  R.  Gupta,  E.  Vredenbregt.  M.  Doerv,  and  other  colleagues 
for  valuable  discussions.  Computing  time  was  provided  by  the  Cornell  National 
Supercomputing  Facility,  which  is  funded  by  NSF  and  IBM. 
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FIGURES 


FIG.  1.  Time  evolution  of  lin  J.  lin  laser  cooling  with  Rb  atoms  (r  =  26.5  ns).  (Top) 
Detuning  =  -  T,  time  interval  between  traces  =  10  /is.  (Bottom)  Detuning  =  -5  T,  time 
interval  =  3  /is. 


FIG.  2.  Average  kinetic  energies  after  lin  J.  lin  cooling  of  Rb  atoms  for  30  /is.  (1130 
r).  See  text  for  definition  of  KEc,  EE  a,  and  Uq. 


FIG.  3.  Calculated  free  particle  velocity  distributions  with  MILC  for  atoms  with  recoil 
energy  3.85  kHz,  r  =  26.5  ns  (as  for  Rb)  for  interaction  time  intervals  of  5  /is  (solid  lines) 
and  50  /is  (dashed  lines)  for  a)  B  =  5  /iT,  S  =  0.42,  6  =  -I\  b)  B  =  5  /iT,  S  =  6.67,  6  - 
-20  T,  c)  B  =  2  /iT,  S  =  0.083,  and  6  =  -I\  and  d)  B  =  7.5  /iT.  S=6.86,  6  =  -20  T.  The 
initial  velocity  distribution  is  flat  to  ±50Vr. 


FIG.  4.  Calculated  populations  of  the  eigenstates  of  the  periodic  light  shift  potential 
for  cases  shown  in  Fig.  2a  (top)  and  2b  (bottom).  The  well  depth  Uo/h  =  1  MHz  m  each 
case.  The  time  intervals  are  0.3  /is  (short  dashed  lines),  5  /is  (solid  lines),  and  50  /is  (long 
dashes).  All  levels  initially  have  a  population  of  1  on  this  scale,  and  subsequently  m?  = 
1/2  levels  become  more  populated,  while  m/r  =  -1/2  levels  fall  below  1.  Band  energies  are 
displayed  at  the  top.  The  distribution  is  non-Maxwellian  for  an  extended  time. 


FIG.  5.  B  field  mixing  coefficients  and  Am/r  =  1  optical  pumping  rates  between 


quantum  levels  in  the  periodic  light  shift  potential,  for  an  F  =  1/2  -*  F  »  3/2  a+  optical 
transition.  Energy  levels  below  the  potential  maxima  (at  1  MHz  in  this  case)  and  bands  are 
shown  for  m p  =  -1/2  (left  in  each  transition  pair)  and  for  mp  =  1/2  (to  the  right  in  each 
pair).  The  width  of  each  arrow  is  proportional  to  the  magnitude  of  the  mixing  coefficient 
or  the  optical  pumping  rate. 
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ABSTRACT 


Sub-Doppler  laaer  cooling  requirea  optical  punping  anong  different¬ 
ly  light  ahifted  ground  atate  sublevels.  In  thia  paper  ve  describe 
a  study  of  the  simplest  possible  angular  momentum  configuration 
that  allow  all  sub-Doppler  cooling  phenomena.  The  Jg  »  1  *  Je  «  o 
angular  momentum  configuration  shows  recoil-limited  cooling  in  the 
two  most  well-known  types  of  polarization  gradients,  magnetically 
induced  laser  cooling,  velocity  selective  resonances,  transient 
cooling,  and  velocity  selective  population  trapping. 


PACS  Numbers:  32.S0.Pj,  42.S0.Vk 


I.  Introduction 


In  all  caaaa  of  sub-Doppler  laser  cooling  (SDLC)  the  role  of 
multiple  ground  state  sublevels,  including  coherences  between 
then  established  by  either  Raman  or  Keenan  processes,  is  of  para¬ 
mount  importance.  The  next  natural  question  is  vhat  is  the  sim¬ 
plest  possible  atomic  level  system  that  would  permit  SDLC.  Since 
multiple  ground  state  sublevels  are  required  for  SDLC,  a  four 
level  0*1  (Jg  «  Je)  transition  will  not  work,  but  J  -  1  *  0  or 
1/2  +  1/2  transitions  will  both  work.  However,  the  1/2  +  1/2 
scheme  will  not  permit  a*  -  ar~  polarization  gradient  cooling 
because  this  requires  atomic  alignment  [1],  so  1  *  0  is  the  sim¬ 
plest  possible  general  scheme  for  SDLC. 

He  have  studied  SDLC  in  the  F  *  1  *  0  transition  of  the  28% 
abundant  87Rb  isotope  in  a  1-D  optical  molasses,  with  and  without 
polarization  gradients,  and  in  weak  and  strong  magnetic  fields 
(strong  means  Zeeman  splitting  larger  than  the  optical  pumping 
rates).  The  single  most  important  characteristic  of  this  transi¬ 
tion  in  light  of  either  linear  or  circular  polarization  is  the 
presence  of  two  unexcitable  ground  state  sublevels  (dark  states) 
as  shown  in  Fig.  1.  Effective  User  cooling  depends  on  very  many 
scattering  events  because  the  momentum  exchange  from  light  scat¬ 
tering  is  typically  much  smaller  than  the  momentum  of  thermal 
atoms.  In  the  usual  J«  J  t  1  transition  schemes,  atoms  are  op¬ 
tically  into  a  cyclic  transition  that  allows  this  required 
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multiple  scattering,  rather  than  into  a  dark  state  that  precludes 
it.  Thus  laser  cooling  on  the  F  «  1  *  0  transition  is  quite 
different  from  the  more  common  experiments  because  transient 
processes  preclude  steady  state  cooling  forces  [2,3]  and  there  is 
no  Doppler  cooling.  On  the  other  hand,  there  is  both  transient 
laser  cooling  and  Velocity  Selective  Population  Trapping  (VSPT) . 

For  the  F  «  1  <*  0  transition  in  a  one-dimensional  field  with 
no  polarization  gradients  there  are  two  dark  states  that  are 
independent  of  atomic  velocity.  However,  with  counterpropagating 
beams  of  different  polarizations  e*  and  C2 *  for  example  circular 
and  linear,  there  is  only  one  velocity-independent  dark  state. 
The  other  two  atomic  ground  state  sublevels  combine  to  form  a 
second  state  that  is  dark  only  at  v  »  0  in  a  semiclassical  pic¬ 
ture,  and  a  third  one  that  is  readily  excited.  This  v  ■  0  dark 
state  in  a  1-D  light  field  with  polarization  gradient  is  a  new 
discovery,  and  will  be  discussed  below.  VSPT  in  this  dark  state 
of  the  1-D  optical  field  can  only  be  eliminated  with  a  strong  B 
field  that  is  not  parallel  to  either  c  for  linear  polarization  or 
to  k  for  circular  polarization.  The  necessity  of  such  e  B  field 
is  easily  understood  by  considering  that  the  e's  have  no  compo¬ 
nents  parallel  to  k,  and  if  k  ware  chosen  as  the  z-axis,  there 
would  be  no  optical  pumping  from  the  Mp  •  0  sublevel  in  the 
absence  of  a  B  field  to  mix  the  states. 
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We  note  three  special  features  of  the  P  ®  1  *  0  transition. 
First,  steady  state  processes  are  easier  to  study  since  non¬ 
steady  state  effects  Bay  be  present  for  a  much  longer  time  for 
larger  Fg  values.  For  example,  optical  pumping  in  circularly 
polarized  light  to  the  steady  state  population  distribution  among 
the  Mp  sublevels  occurs  before  steady  state  cooling  on  a  cyclic 
transition  begins  [3].  Second,  calculations  of  the  damping  force 
or  the  velocity  distribution  are  simplified  because  of  the  rela¬ 
tively  small  Hamiltonian  matrix.  Third,  Fe  •  0  requires  that  the 
ground  states  can  be  optically  coupled  to  only  one  excited  state 
so  that  ground  state  coherences  are  not  diluted  by  multiple 
excited  states. 
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In  section  II  below  we  present  a  detailed  discussion  of 
several  special  phenomena  of  laser  cooling  in  the  F  ■*  1  *  0  tran¬ 
sition.  In  section  III  we  describe  our  apparatus  and  present  the 
experimental  results,  including  observation  of  VSPT.  In  section 
IV  we  summarize  our  findings. 


II.  Sub-Doppler  Cooling  on  in  F  ■  1  *  0  Transition. 

The  special  characteristics  of  laser  cooling  on  an  F  *  1  #  0 
transition  can  be  used  to  elucidate  many  features  of  SDLC.  We 
first  consider  laser  cooling  in  either  the  -  a r  or  the  lin  i 
lin  configurations  (4]  with  a  strong  D  field  applied  perpendicu- 
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lar  to  the  optical  k  vectors  (for  the  lin  i  lin  case,  B  must  not 
be  parallel  to  either  of  the  polarizations  c).  In  these  cases 
each  traveling  wave  laser  beam  can  induce  both  n  and  <r  transi¬ 
tions  in  the  frame  with  quantization  axis  chosen  along  the  strong 
magnetic  field.  Then  both  polarization  schemes  permit  an  energy- 
exchange-based  description  of  laser  cooling  toward  v  *  0 

(Sisyphus  cooling) . 

To  understand  this  description  we  view  the  optical  field  in 
either  polarization  scheme  (with  B  at  ±45°  to  the  c's  of  lin 
lin)  as  two  orthogonal,  linearly  polarized  standing  waves  with 
nodes  spatially  displaced  from  one  another  as  shown  in  Fig.  2a. 
Such  a  field  is  symmetric  for  rotation  about  k,  so  we  choose  one 
of  the  c's  along  B.  Thus  one  standing  wave  field  induces  n  tran¬ 
sitions  and  the  other  one  induces  <r  transitions. 

We  consider  light  of  sufficiently  low  intensity  and  large 
detuning  so  that  the  excited  state  population  is  small.  Atoms 
are  then  subject  to  the  conservative  force  from  the  light  shift 
potential  (calculated  by  adiabatically  eliminating  the  excited 
state).  We  consider  the  conservative  motion  of  the  atoms  in  the 
sinusoidal  light  shift  potentials,  coupled  with  the  optical  pump¬ 
ing  that  provides  the  required  irreversibility  to  produce  damping 
forces  on  the  atomic  motion,  as  shown  schematically  in  Fig.  2b. 
Mote  that  tarn  strengths  of  the  <r  transitions  are  smaller  than 
those  of  u  transitions  for  the  same  linear  polarisation  amplitude 
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|c|,  so  the  potential  hills  are  correspondingly  smaller.  Since 
the  potential  hills  associated  with  the  light  shift  for  the  stan¬ 
ding  wave  of  one  polarization  are  displaced  from  those  of  the 
other,  optical  pumping  can  switch  the  atomic  populations  among 
the  ground  state  sublevels  so  that  moving  atoms  always  climb  the 
potential  hills  [5]. 


This  is  indeed  what  happens  when  the  detuning  from  atomic 
resonance,  5  »  «iaeer  “  watom  >  °*  Atoms  that  enter  the  light 
field  near  antinode  of  standing  wave  are  quickly  optically 
pumped  out  of  the  sublevel  excited  by  c*.  Such  atoms  can  not  be 
excited  by  C2  because  they  are  at  the  node  for  this  polarization 
where  the  excitation  rate  is  zero  (as  is  the  light  shift) .  Their 
transverse  velocity  carries  them  towards  the  antinode  of  C2,  and  • 

in  the  process  they  increase  their  internal  energy  because  of 
positive  light  shift,  and  correspondingly  decrease  their  kinetic 
energy  and  slow  down.  Near  the  antinode  of  C2  where  the  light 
intensity  and  light  shift  are  largest,  the  slowed  atoms  are  opti¬ 
cally  pumped  into  the  state  excited  by  e^,  but  they  are  now  at 
the  node  of  this  standing  wave.  The  energy  loss  process  repeats, 
and  the  atomic  sample  is  cooled. 


In  contrast  to  the  energy  based  Sisyphus  picture  described 
above,  laser  cooling  in  these  polarisation  configurations  can 
also  be  viewed  in  the  momentum-based,  Velocity  Selective  Reso¬ 
nance  (VSR)  picture.  In  this  case  we  necessarily  choose  a  des- 
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cription  of  the  light  field  as  counterpropagating  travelling 
waves  instead  of  standing  waves.  When  B  is  strong  enough  for  the 
Zeeman  splitting  to  be  larger  than  the  optical  pumping  rate  rp, 
and  thereby  make  B  a  suitable  choice  for  the  quantization  axis, 
each  light  beam  (not  each  standing  wave)  in  either  the  a+  -  a~  or 
the  lin  x  lin  configuration  can  induce  both  n  and  <r  transitions. 
Thus  there  can  be  VSR,  comprising  excitation  by  one  beam  followed 
by  stimulated  emission  by  the  other,  back  to  the  original  ground 
state  sublevel  as  shown  by  the  dashed  arrows  in  Fig.  3.  Such  a 
sequence  can  produce  strong  velocity  damping  toward  the  resonance 
velocity,  which  is  v  *  0  in  this  case,  because  the  initial  and 
final  states  have  the  same  energy  (same  state)  [6,7]. 

This  VSR  picture  provides  insight  beyond  the  Sisyphus  pic¬ 
ture  because  there  are  also  VSR's  between  pairs  of  sublevels  of 
different  energy.  These  are  levels  that  are  Zeeman  split  by  wz  ■ 
9F^Bb  [6,7] ,  and  Raman  transitions  between  them  are  shown  by  the 
solid  arrows  in  Fig.  3.  Such  Raman  transitions  are  resonant  when 
2k-  v  »  ±u>z,  and  are  resolvable  when  7p  <  wz*  Our  data  presented 
below  show  strong  sub-Doppler  cooling  to  both  v  -  0  and  v  ■ 

±tdz/2k. 


III. A.  Experimental  Apparatus 


Most  of  our  experimental  setup  has  been  described  in  previ¬ 
ous  work  [6,8]  and  is  only  briefly  described  here.  Rb  atoms 
emerge  from  a  150  °c  oven  with  a  horizontal  slit  aperture  0.06  mm 
high  by  2  mm  wide  (see  Fig.  4) .  A  horizontal  atomic  beam  is 
formed  by  a  vertical  slit  2  mm  high  by  0.06  mm  wide,  35  cm  away 
from  the  oven.  The  atomic  beam  is  crossed  at  90°  by  a  pair  of 
counterpropagating ,  horizontal  laser  beams  just  beyond  the  verti¬ 
cal  slit,  and  the  polarization  of  the  retroreflected  beam  can  be 
modified  by  a  quarter  wave  plate  in  front  of  the  mirror.  Three 
pairs  of  square  Helmholtz  coils  are  used  for  controlling  the  B 
field  in  the  interaction  region.  The  atomic  beam  profile  paral¬ 
lel  to  the  direction  of  the  optical  k  vectors  is  measured  with  a 
scanning  hot  platinum-tungsten  wire,  25  um  in  diameter,  1.3  m 
away  from  the  interaction  region. 

A  20  mW  Sharp  model  LT024  diode  laser  is  side-locked  to  a 
Doppler-broadened  signal  from  a  Rb  cell  at  room  temperature. 

Light  from  this  laser  is  used  to  optically  pump  the  87Rb  atoms 

\ 

from  the  5S1/2(F»2)  to  the  5S1/2(F»1)  ground  hfs  state  before 
they  reach  the  interaction  region. 

Light  from  a  second  diode  laser,  a  35  mW  Sharp  model  LT025, 
is  split  into  two  beams,  and  the  weaker  one  passes  through  an 
80  MHz  AOM  and  then  to  another  Rb  vapor  cell.  A  cross-over  reso- 


nance  in  its  saturated  absorption  signal  is  used  to  lock  the 
laser  near  the  F  *  1  *  0  transition  in  87Rb  near  X  *  780  nm,  and 
facilitates  tuning  on  either  side  of  the  atomic  resonance.  About 
10**5  of  the  laser  light  is  fed  back  to  the  laser  by  a  blazed 
grating  (see  Fig.  4)  to  reduce  the  spectral  width  to  ~1  MHz  [9]. 
The  laser  beam  has  its  spatial  intensity  profile  flattened  to  a 
few  percent  by  a  tilted  etalon  [10],  and  is  then  expanded  to  23 
mm  wide  by  5  mm  high.  A  mechanical  shutter  interrupts  this  beam, 
and  many  systematic  effects  have  been  eliminated  by  subtracting 
data  from  measurements  with  the  laser  beam  off  from  those  with 
the  beam  on. 
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Figure  5  shows  the  measured  atomic  beam  profiles  with  B  per¬ 
pendicular  to  the  k  vectors  for  (a) ,  the  case  of  cr+  -  <r~  and  (b) 
lin  i  lin  with  B  at  45°  to  the  polarization  vectors.  In  each 
case  B  is  perpendicular  to  the  horizontal  atomic  and  laser  beams. 

The  top  trace  in  each  figure  shows  strong  sub-Doppler  cooling  to 

\ 

v  »  0,  and  a  velocity  width  of  only  -  5.5  cm/s  (Doppler  limit  * 
13  cm/s  in  1-D  optical  molasses) .  Both  the  top  traces  were  taken 
with  large  B  (B  >  1  G) ,  and  both  also  show  the  VSR  peaks  at  v  * 
0.  As  B  is  reduced  (lower  traces),  these  VSR  peaks  move  to  lower 
velocities,  and  Fig.  6  shows  the  separation  of  these  peaks  vs. 

| B |  for  each  case. 
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In  the  semiclassical  approximation  the  "velocity  capture 
range"  of  this  cooling  process  is  determined  by  the  lifetime  of 
the  ground  state  populations  and/or  coherences  produced  by  the 
VSR's.  Atoms  are  optically  excited  out  of  these  coherences  at  a 
rate  rp  -  fp  cos2(kx)  where  fp«  2s(y/2) / (3L  +  s ),  x  ■  1/7  is  the 
excited  state  lifetime  (27  ns  for  this  Rb  transition) ,  L  ■  1  + 
(26/7) 2 /  and  the  saturation  parameter  for  a  single  light  beam  is 
s  «  I/ISat»  where  I8at  *  nhc/3X3x  (s  1.6  mW/cm 2  for  the  X  *  780 
nm  transition  in  Rb) .  The  factor  of  3  in  the  denominator  of  frp 
comes  from  the  Clebsch-Gordan  coefficient  for  the  1  *  0  transi¬ 
tion  for  each  value  of  Mp,  and  the  factor  of  2  in  the  numerator 
comes  from  the  superposition  of  two  orthogonal  light  beams.  For 
Fig.  5a,  s  #  1.33  and  8  «  2n  x  12  MHz  whence  fp/2 n  «  150  kHz  and 
the  average  value  over  a  wavelength  ?p  is  half  of  that. 


I 


I 


(4r) 


» 


> 


> 


>  • 


Under  optimum  cooling  conditions,  the  expected  width  of  the 
velocity  distribution  is  -  ?p/k  *  6  cm/s,  about  equal  to  the  va-  ) 

lue  observed  for  the  v  *  0  peak.  The  widths  of  the  v  *  0  peaks, 

Av,  are  limited  by  the  width  of  the  longitudinal  velocity  distri¬ 
bution  Av{  through  Av/v  -  (Av^/v/)  -  1/2  for  a  thermal  beam,  and  i 

their  widths  are  indeed  about  half  their  distances  from  the  v  »  o 

peak.  They  begin  to  merge  into  it  when  their  velocity  is  -  10 
cm/s,  corresponding  to  |B|  a  0.3  G  for  gp  ■  1/2.  similar  conclu-  > 

sions  can  be  drawn  from  the  data  of  Fig.  5b.  Furthermore,  at  low 
B  fields  some  states  are  only  weakly  excited  and  transient 
cooling  effects  begin  to  appear  [2].  i 
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The  two  bottom  traces  in  Fig.  5,  where  |Bj  *  0,  reveal  the 
velocity  selective  population  trapping  (VSPT)  discussed  above. 
Atoms  travelling  in  either  light  field  of  Fig.  5  with  nearly  zero 
transverse  velocities  spend  all  their  time  in  light  of  one  parti¬ 
cular  polarization  as  they  ride  along  a  particular  phase  of  the 
standing  wave.  For  example,  in  lin  J-  lin  this  may  be  circular 
along  one  path  across  the  standing  wave,  or  linear  along  another 
path  that  is  X/8  away.  In  the  absence  of  a  B  field  the  quantiza¬ 
tion  axis  is  chosen  along  the  local  polarization.  Atoms  are  then 
optically  pumped  to  a  dark  state  of  this  basis  (Fig.  1)  and  re¬ 
main  uncoupled  to  the  excited  state  by  the  light.  These  slow 
atoms  are  subject  to  VSPT  and  thus  do  not  experience  momentum 
diffusion.  They  could  therefore  display  an  arbitrarily  narrow 
velocity  distribution  with  a  long  enough  interaction  time.  (Note 
that  this  is  not  the  same  as  the  "coherent"  population  trapping 
of  Ref.  11  which  depends  on  the  quantized  center-of-mass  motion 
of  the  atoms.)  Atoms  are  cooled  to  the  low  velocities  appropri¬ 
ate  for  VSPT  by  transient  polarization  gradient  cooling,  and 
these  form  the  sub-Doppler  peaks  at  the  bottom  of  Fig.  5  with 
B  «  0.  When  atoms  traverse  either  light  field  of  Fig.  5  with  a 
higher  velocity,  they  experience  a  constantly  changing  polariza¬ 
tion,  and  could  thus  be  pumped  out  of  the  velocity  selective 
"dark"  states  of  Fig.  1. 
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We  have  tested  this  hypothesis  of  VSPT  at  |§|  -  o  several 
ways.  First,  we  found  that  the  cold  atom  peaks  occur  over  a  wide 
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range  of  detuning,  both  blue  and  red,  although  they  are  stronger 
in  the  blue  case  because  there  is  residual  polarization  gradient 
cooling  to  feed  the  dark  state.  Second,  we  used  an  aperture  on 
the  laser  bean  to  shorten  the  interaction  tine,  and  found  that 
the  cold  aton  peaks  require  a  few  optical  punping  tines  1/fp  to 
appear.  Third,  we  note  that  fp  m  y/80  for  the  conditions  of  Fig. 
5,  so  that  optical  punping  occurs  in  about  2.5  4s.  VSPT  atons 
that  nove  less  than  1/4  in  this  tine  will  have  to  wait  longer 
than  this  tine  to  be  punped  out  of  the  dark  state  because  they 
will  not  be  subject  to  a  different  polarization.  Thus  there  is 
also  a  cold  peak  even  in  the  absence  of  a  B  field  with  a  velocity 
distribution  of  width  -  9  cn/s,  about  50%  larger  than  the  large  B 
field  peaks,  just  as  shown  at  the  botton  of  Fig.  5. 
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In  addition  to  all  the  experiments  described  above  with  the 
two  types  of  polarization  gradients,  we  have  also  studied  SDLC  in 
this  F  »  1  *  0  transition  in  a  standing  wave  of  purely  circularly 
polarized  light.  The  presence  of  a  weak  transverse  B  field  then 
produces  nagnetically  induced  laser  cooling  [8].  At  low  B  field 
we  observed  a  very  narrow  peak  near  v  »  0,  and  at  higher  fields 
we  saw  this  split  into  two  peaks  with  sub-Doppler'  widths  centered 
at  vr  ■  ±0)2/ 2k  just  as  in  previous  experinents  using  sore  compli¬ 
cated  transitions  [6].  We  have  also  observed  higher  order  VSR's 
where  the  resonance  velocity  satisfies  vr'  •  ±wz/4k,  correspon¬ 
ding  to  redistribution  of  two  photons  between  the  two  laser  beans 
that  torn  the  optical  nolasses. 


» 


IV.  Conclusions. 


We  have  performed  several  other  experiments  on  this  simplest 
possible  SDLC  scheme  to  assure  that  it  displays  all  the  features 
%re  expect.  For  example,  in  the  lin  J-  lin  case  with  B  along  one 
of  the  polarization  vectors,  ve  find  that  both  the  cold  peaks  of 
atoms  indeed  move  at  the  resonance  velocity  vr  ■  tog /2k  [12]. 

Thus  we  conclude  that  the  F  *  1  «*  0  transition  is  the  para¬ 
digm  of  SDLC,  displaying  many  of  its  special  features  in  addition 
to  uK  temperatures,  and  holding  claim  to  being  the  simplest  pos¬ 
sible  angular  momentum  configuration  with  this  capability.  In 
the  future  we  plan  more  extensive  studies  of  VSPT,  the  higher 
order  resonances,  and  other  aspects  of  this  extraordinary  laser 
cooling  system. 
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FIGURE  CAPTIONS 


Fig.  l.  Ths  energy  level  scheme  for  an  F  ■  1  *  0  transition 
showing  allowed  excitation  for  blue  tuned  <r+  light  (a)  and  n 
light  (b).  In  (a)  the  Mp  -  0  and  1  states  are  dark,  whereas  in 
(b)  the  Mp  -  i  1  states  are  dark.  Linearly  polarized  light  per¬ 
pendicular  to  z  would  excite  both  the  Hp  ■  4-1  and  -1  sublevels 
leaving  Mp  -  0  as  one  dark  state,  while  an  uncoupled  superposi¬ 
tion  of  the  Mp  «  ±1  states  would  be  the  other. 

Fig.  2a.  The  standing  wave  fields  of  either  a*  -  or  lin 
-L  lin  polarization  configurations.  Both  optical  fields  have  the 
sane  xr/2  spatial  phase  shift.  The  only  difference  between  the 
the  two  polarization  schemes  is  the  relative  temporal  phase  of 
the  linearly  polarized  standing  waves.  This  phase  difference  is 
zero  for  <r+  -  <r~  and  ±n/2  (i.e.,  time  lag  of  *K/2<*>iaser)  for  lin 
-1  lin.  The  selection  rules  are  the  same.  If  B  is  along  one  of 
the  c's  in  the  lin  x  lin  case,  one  standing  wave  is  parallel  to 
the  quantization  axis  along  B  so  it  induces  n  transitions,  and 
the  other  one  is  perpendicular  to  B  so  it  induces  both  o+  and  a’ 
transitions.  In  the  cr+  -  case,  choosing  i  along  B  yields  the 
seme  selection  rules  as  the  lin  x  lin  case  (see  Ref.  7).  2b 
shows  the  spatial  dependence  of  the  light  shift  for  the  field  of 
as  when  the  quantitation  axis  is  chosen  perpendicular  to  &  (at 
*45°  to  the  c  fields  for  lin  ±  lin) .  Since  the  strength  of  the  n 
transition  is  larger,  the  light  shift  is  larger  for  the  same  |c  j . 


Fig.  3.  Raman  transitions  among  the  Mp  sublevels  of  the  F  * 
1  +  0  transition  in  a  magnetic  field.  VSR  at  v  ■  0  are  induced 
in  atoms  that  return  to  their  original  sublevel  (dashed  arrows) , 
and  at  finite  velocity  vr  ■  wj/2k  for  atoms  that  are  transferred 
between  states  of  different  Mp  (solid  arrows) . 

Fig.  4.  Schematic  diagram  of  the  apparatus. 

Fig.  5.  Measured  atomic  beam  profiles  for  5/2 n  -  +12  MHz  (e 
2 r)  for  the  case  of  (a)  with  s  *  1.33,  and  (b)  lin  «i  lin 
with  s  *  0.66.  The  highest  B  values  are  at  the  top,  decreasing 
to  B  *  0  at  the  bottom.  Note  that  the  VSR  signals  on  each  side 
of  the  central  peak  are  broadened  by  the  longitudinal  velocity 
distribution. 

Fig.  6.  Velocity  of  the  VSR  peaks  in  Fig.  5  vs  B.  The  so¬ 
lid  dots  are  for  <r+  (Fig.  5a)  and  the  open  squares  are  for  lin  -1 
lin  (Fig.  5b).  The  straight  lines  represent  the  resonance  case, 
v  «  ug/2k  appropriate  in  each  case  because  the  light  induces  both 
it  and  a  transitions  in  the  strong  transverse  B  field. 
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ABSTRACT 

We  report  a  new  type  of  sub-Doppler  laser  cooling  in  1-D 
using  neither  polarization  gradients  nor  magnetic  fields.  Two 
laser  frequencies  couple  each  of  the  two  ground  state  hfs  sub- 
levels  to  a  common  excited  state  of  85Rb,  thus  forming  a  Asystem. 
We  observe  sub-Doppler  cooling  when  the  light  is  tuned  blue  of 
resonance,  efficient  heating  with  red  detuning,  deflection  from  a 
rectified  dipole  force,  and  velocity  selective  resonances  associ¬ 
ated  with  Raman  transitions.  A  simplified  semiclassical  calcula¬ 
tion  agrees  qualitatively  with  our  measurements. 


PACS  numbers  32.80.Pj,  42.50.Vk 
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All  sub-Doppler  laser  cooling  (SDLC)  methods  can  be  charac¬ 
terized  by  a  damping  force  derived  from  a  spatially  varying  light 
field  (e.g.,  from  superposition  of  two  counterpropagating  laser 
beams)  [1,2].  In  all  cases  studied  to  date,  motion  of  the  atoms 
prevents  their  internal  state  distribution  from  reaching  steady 
state.  Deviation  from  the  steady  state  is  countered  by  the  ir¬ 
reversible  process  of  optical  pumping  (spontaneous  emission)  that 
tends  to  restore  the  steady  state  distribution.  The  time  lag  of 
this  process  leads  to  non-adiabatic  response  of  moving  atoms  to 
the  standing  wave  field,  and  under  suitable  conditions,  a  velo¬ 
city  dependent  damping  force.  In  the  past  polarization  gradients 
[1,2]  or  magnetic  fields  [3]  were  used  to  create  the  spatial 
dependence  of  the  atomic  steady  state.  In  this  letter  we  report 
a  new  SDLC  scheme  that  requires  neither  polarization  gradients 
nor  magnetic  fields. 

Our  new  scheme  employs  two  standing  waves  with  the  same  lin¬ 
ear  polarization  but  different  laser  frequencies,  Wi  and  u2,  cor¬ 
responding  to  excitation  of  the  two  ground  state  hyperfine  levels 
of  85Rb  (separated  by  w^fS  a  2n  x  3.036  GHz).  When  u>i,  u2,  and 
6  *  (w2  -^i)  -  Whfs  are  properly  chosen,  the  time  lag  described 
above  can  also  damp  the  motion  of  moving  atoms  because  of  the 
non-adiabatic  response  of  the  populations  of  levels  |l>  and  |2>. 
Recently  there  has  been  considerable  theoretical  interest  in  the 
optical  forces  on  three  level  atoms  with  such  a  A  configuration 
of  levels  [4-12]  but  only  a  few  experiments  [13-15]. 

SDLC  is  often  described  in  a  picture  called  Sisyphus  cooling 
(because  of  the  close  analogy  with  a  Greek  myth  [16])  where  atoms 
lose  kinetic  energy  when  they  move  in  the  potentials  caused  by 
the  light  shifts  of  the  laser  field.  Such  descriptions  exist  at 
low  intensity  for  the  steady  state  processes  of  polarization 
gradient  [1,16,17]  and  magnetically  induced  cooling  [3],  as  well 
as  transient  laser  cooling  processes  [18].  With  two  frequencies 
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present,  a  similar  energy  loss  or  Sisyphus  picture  is  appropriate 
when  the  relative  spatial  phase  of  the  standing  waves  is  consi¬ 
dered. 

Consider  the  case  where  the  intensity  maxima  of  the  standing 
wave  at  u>i  correspond  to  the  minima  of  u>2  as  shown  in  Fig.  la, 
and  both  wj/s  are  larger  than  the  corresponding  atomic  resonance 
frequencies  as  shown  in  Fig.  lb.  The  light  shift  of  each  ground 
state  sublevel  |i>  is  dominated  by  the  nearly  resonant  frequency 
and  hence  is  positive.  Atoms  moving  through  the  node  of  <*>! 
are  optically  pumped  to  1 1>  and  must  increase  their  potential 
energy  to  move  further  as  shown  in  Fig.  la.  This  comes  at  the 
expense  of  their  kinetic  energy,  and  they  slow  down  as  they 
approach  the  top  of  the  potential  hill  at  the  antinode  of 
However,  they  are  most  likely  to  be  optically  pumped  to  |2>  at  a 
point  where  the  intensity  of  light  at  is  high,  and  then  they 
begin  climbing  another  hill.  Repetition  of  this  process  rapidly 
converts  kinetic  energy  into  potential  energy  which  is  then  dis¬ 
sipated  into  the  fluorescent  light  of  the  spontaneous  decays,  -and 
atoms  are  cooled.  This  process  obviously  produces  heating  if  the 
two  frequencies  are  red  of  resonance  thereby  producing  negative 
light  shifts,  and  its  effectiveness  is  compromised  when  the  stan¬ 
ding  wave  phase  difference  is  different  from  n. 

Most  of  our  experimental  setup  has  been  described  in  previ¬ 
ous  work  [3,17,18]  and  is  only  briefly  described  here.  Rb  atoms 
emerge  from  a  150  °C  oven  with  a  horizontal  slit  aperture  0.06  mm 
high  by  2  mm  wide.  A  horizontal  atomic  beam  is  formed  by  a  ver¬ 
tical  slit  2  mm  high  by  0.06  mm  wide,  35  cm  away  from  the  oven. 
The  atomic  beam  is  crossed  at  90°  just  after  the  second  slit  by  a 
laser  beam  containing  two  frequencies  that  is  retroreflected  by  a 
movable  mirror  to  form  two  horizontal  standing  waves  with  the 
same  polarization  but  different  frequencies.  The  transit  time  of 
atoms  with  the  average  velocity  350  m/s  in  the  23  mm  long  by  5  mm 
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high  laser  beam  is  about  65  ms.  The  Earth's  magnetic  field  is 
carefully  cancelled  to  less  than  10  mG  by  square  Helmholtz  coils. 
We  measure  the  atomic  beam  profile  parallel  to  the  direction  of 
the  optical  k  vectors  with  a  scanning  hot  platinum-tungsten  wire, 
25  m  in  diameter,  1.3  m  away  from  the  interaction  region. 

Linearly  polarized  light  from  a  35  mW  Sharp  model  LT025 
diode  laser  is  split  into  two  beams,  and  the  weaker  one  passes 
through  a  80  MHz  AOM  and  then  to  a  Rb  vapor  cell.  A  crossover 
resonance  in  the  saturated  absorption  signal  is  used  to  lock  the 
laser  near  the  F  =  3  ■»  3  transition  in  85Rb  near  X  =  780  nm,  and 
facilitates  tuning  on  either  side  of  the  atomic  resonance.  About 
10~5  of  the  stronger  beam  is  fed  back  to  the  laser  by  a  blazed 
grating  to  reduce  the  spectral  width  to  -l  MHz.  The  two  frequen¬ 
cies  are  generated  by  passing  the  light  through  an  EOM  [19].  A 
mechanical  shutter  interrupts  this  beam,  and  many  systematic 
effects  have  been  eliminated  by  subtracting  laser  beam  off  from 
laser  on  measurements. 

The  spatial  phase  difference  between  the  two  light  shift 
potentials  is  of  paramount  importance  in  this  experiment.  The 
standing  waves  of  the  two  frequencies  have  nodes  at  the  surface 
of  the  retroref lecting  mirror,  but  their  relative  phase  *  = 
2z(u>i  -  c>2)/c  varies  with  distance  z  from  the  mirror  because  of 
their  frequency  difference;  for  <«>i  -  u>2  »  2n  x  3  GHz  the  spatial 
period  of  this  variation  is  a  5  cm.  Since  this  is  much  larger 
than  the  0.06  mm  atomic  beam  diameter,  atoms  see  a  fixed  value  of 
the  phase  between  the  standing  waves.  This  can  be  varied  by 
translation  of  the  retroref lecting  mirror. 

Figure  2a  shows  the  results  of  our  measurements  of  SDLC 
using  two  optical  frequency  standing  waves  for  the  case  of  •  ■  n 
(optimum  cooling)  and  •  -  0  obtained  by  moving  the  retroreflec- 
ting  mirror  about  2.5  cm.  Even  at  t  »  0  there  is  a  significant 
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signal  of  cooled  atoms,  but  this  arises  from  two  additional  cool¬ 
ing  mechanisms  operating  in  this  experiment.  First,  the  atomic 
coherences  produced  by  stimulated  Raman  transitions  between  the 
states  |l>  and  |2>  for  5  *  0,  neglected  in  this  Sisyphus  descrip¬ 
tion,  produce  velocity  selective  resonances  (VSR)  at  v  «  0  as 
described  below.  Second,  for  the  F  *  3  3  transition  driven  by 
linearly  polarized  light,  the  excitation  of  the  Mp  *  0  sublevel 
is  forbidden,  and  transient  cooling  into  this  dark  state  contri¬ 
butes  to  this  residual  cooling  [18]. 

We  have  tested  that  this  residual  peak  of  cold  atoms  is 
mostly  the  result  of  transient  cooling  by  varying  the  interaction 
time  with  an  aperture  on  the  laser  beams.  We  found  that  it  is 
produced  in  a  time  comparable  to  the  optical  pumping  time  l/yp. 
For  this  test  we  used  a  large  enough  value  of  £  so  that  the  VSR 
peaks  are  well  resolved  from  the  v  *  0  signal,  (see  below),  and 
the  remaining  signal  at  *  =  0  is  only  from  transient  cooling. 
Figure  2b  shows  the  4  dependence  of  the  height  of  the  cooled  atom 
peak  at  v  =  0  with  large  6. 

The  Sisyphus  description  of  these  experiments  given  above  is 
limited  because  it  is  based  only  on  the  populations  of  states, 
and  neglects  coherences  that  are  associated  with  Raman  processes 
between  the  two  long  lived  ground  states.  It  is  well  known  that 
these  coherences  have  important  effects  on  the  damping  forces, 
and  in  previous  experiments  we  have  shown  that  damping  to  v  «  0 
[20]  and  to  finite  velocities  [20,21]  depends  on  Raman  resonances 
between  Zeeman  shifted  states  that  redistribute  light  between 
counterpropagating  beams.  These  VSR  are  always  present  in  all 
known  cases  of  steady  state  SDLC  [22,23]. 

The  direct  analog  of  the  magnetic  VSR  can  be  seen  with  two 
laser  frequencies  when  £  •  (c*>2  -  Uj.)  -  u^fs  *  0  [20].  Then  the 
condition  for  VSR  is  satisfied  at  the  resonance  velocity  vr  - 
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±  6 /2k.  Figure  3a  shows  the  measured  atomic  beam  profiles  for 
different  values  of  5  and  Fig.  3b  shows  the  velocity  associated 
with  the  side  peaks  of  Fig.  3a  vs  6.  Two  side  peaks  appear  for 
each  value  of  5  because  atoms  at  the  resonance  velocity  vr  can 
absorb  either  or  u>2  from  either  direction  of  the  laser  beam, 
depending  on  the  sign  of  their  velocity.  Of  course,  stimulated 
emission  into  the  other  beam  completes  the  Raman  transition. 

The  widths  of  the  v  #  0  peaks,  Av,  are  limited  by  the  width 
of  the  longitudinal  velocity  distribution  Av*  through  Av/v  - 
(Av*/v*)  -  1/2  for  a  thermal  beam,  and  their  widths  are  indeed 
about  half  their  distances  from  the  v  =  0  peak.  They  begin  to 
merge  into  it  when  their  velocity  is  -  20  cm/s,  corresponding  to 
6  a  0.25  MHz.  It  must  be  emphasized  that  the  widths  of  these 
peaks  are  all  below  the  Doppler  limit  (10  cm/s):  atoms  undergo 
SDLC  to  the  velocity  vr. 

Up  to  this  point  we  have  considered  only  the  special  cases  * 
*  0  or  n ,  and  the  same  sign  for  the  detuning  *  <*>i  -  <>>atom(i) 
of  the  two  individual  laser  frequencies  (see  Fig.  1).  We  now 
consider  cases  where  these  conditions  are  relaxed.  In  addition 
to  SDLC,  we  show  that  various  values  of  6*  and  52  can  produce 
qualitatively  new  effects.  In  particular,  when  6*  and  62  have 
opposite  signs,  the  light  shifts  for  the  two  ground  states  ji> 
and  1 2>  are  opposite,  for  f  »  0  or  it,  partial  cancellation  of 
the  kinetic  energy  change  when  the  atoms  move  on  these  oppositely 
light  shifted  potentials  severely  reduces  the  heating  or  cooling 
forces,  but  for  other  values  of  t  this  cancellation  is  less  com* 
plate,  and  there  remains  a  "rectified**  dipole  force  [11,13,15]. 

Figure  4  shows  schematically  how  this  works  if  we  choose  for 
simplicity  equal  probability  for  the  optical  pumping  processes 
|1>  «*  1 2>  and  }  2>  <♦  \  1>.  Atoms  ere  then  optically  pumped  to  the 
state  with  the  lowest  standing  wave  intensity  as  shown  by  the 
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arrows  in  Fig.  4.  For  velocity  v  snail  enough  so  that  the  aver* 
age  optical  punping  tine  l/7p  is  less  than  Az/v,  where  Az  is  the 
standing  wave  displacement,  then  the  optical  force  (slope  of  the 
potential  curve  to  which  atons  are  optically  pumped)  has  a  fixed 
direction  independent  of  velocity  when  averaged  over  a  wave¬ 
length.  We  therefore  expect  a  maximum  force  near  v  ■  0,  with  the 
sane  sign  for  v  <  0  or  v  >  0.  The  force  should  decrease  with  a 
velocity  width  -  ?pAz  a  ?p/k. 

In  our  experiment  the  atoms  emerge  from  the  second  slit  in  a 
velocity  range  of  -  ±1  m/s,  much  larger  than  the  Mcapture  range" 
7p/k  of  -  10  cm/s.  Atoms  with  a  large  transverse  velocity  are 
thus  unaffected  by  the  deflection  force,  while  slower  atoms 
undergo  a  deflection  to  the  right  or  left  depending  on  the  sign 
of  Az.  Changing  the  sign  of  Az  can  be  accomplished  either  by  ex¬ 
changing  the  signs  of  the  8j/s  or  by  moving  the  retroref lection 
mirror,  and  we  have  done  both.  The  data  of  Fig.  5a  and  5b 
clearly  show  such  deflections  of  the  atoms  for  opposite  signs  of 
Az,  and  also  provide  a  measure  of  the  effective  velocity  range 
for  this  rectified  force. 

We  have  done  semiclassical  calculations  of  the  force  on  mov¬ 
ing  ato&s,  starting  with  the  optical  Bloch  equations  to  find  the 
density  matrix.  Using  only  the  transfer  of  population  between 
states  1 1>  and  |2>,  and  neglecting  any  coherences  that  may  be 
produced  in  the  optical  excitation,  we  have  calculated  the  force 
for  the  many  experiments  described  above.  For  the  data  of  Fig.  2 
we  find  the  usual  dispersion-shaped  velocity  dependence  with 
strength  and  capture  range  that  correspond  well  with  our  measure- 
manta.  For  the  data  of  Fig's.  5a  and  5b,  we  show  in  Fig's.  Sc 
and  5d  the  anticipated  single  peak  (or  dip)  in  the  force  vs  velo¬ 
city  curve,  with  no  zero  crossing  and  a  width  -  yp/k.  such  model 
calculations  can  not  be  used  for  the  data  of  Fig.  3  because  the 
VSR  arise  from  coherences  that  are  not  included.  We  plan  to 
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extend  the  semiclassical  methods  ve  have  developed,  that  include 
the  multiple  magnetic  sublevels  [22,23],  to  two  frequencies  to 
model  these  experiments. 

Cooling  atoms  with  two  frequencies  is  a  fundamentally  new 
SDLC  scheme  that  requires  neither  polarization  gradients  nor 
magnetic  fields.  Its  damping,  maximum  force,  capture  range,  and 
diffusion  are  comparable  to  those  of  many  previously  studied  SDLC 
schemes  such  as  polarization  gradients  [1,2]  and  magnetically 
induced  laser  cooling  [3].  However  there  are  additional  free 
parameters  that  provide  flexibility  for  further  study  of  SDLC 
processes.  These  are  the  detunings  5^  and  62  >  saturation 
parameters  s^  and  83,  and  the  phase  between  the  standing  waves  4. 
It  also  provides  a  better  insight  into  particular  interpretations 
of  SDLC  such  as  the  Sisyphus  picture.  We  plan  to  extend  our 
study  of  this  system  to  include  polarization  differences  between 
the  two  beams,  polarization  gradients  of  each  beam,  and  the 
effect  of  magnetic  fields. 


I 


(+• 


1 


ft 


> 


»  • 


Supported  by  NSF,  ONR,  AFOSR,  and  CAPES  (Brazil) . 


1 


I 


8 


© 


•  •  • 


REFERENCES 

1.  J.  Dalibard  and  C.  Cohen-Tannoud j i ,  J.  Opt.  Soc.  An.  ££,  2023 
(1989). 

2.  F.  Ungar  et  al.,  J.  Opt.  Soc.  An.  fi£»  2058  (1989). 

3.  B.  Sheehy  et  al.,  Phys.  Rev.  Lett.  £4,  858  (1990). 

4.  N.  Lu  et  al.,  Phys.  Rev.  A35.  5088  (1986) 

5.  J.  Javanainen,  Phys.  Rev.  Lett.  £4,  519  (1990). 

6.  S.  Chang  et  al..  Opt.  Comnun.  12*  19  (1990). 

7.  M.  Prentiss  et  al.,  Opt.  Lett.  1£,  1695  (1991). 

8.  Y.  Rozhdestvenskii  and  N.  Jakobsen,  Sov.  Phys.  JETP  12,  936,(1991). 

9.  D.  Kosachiov  et  al- ,  Opt.  Commun.  ££,  209  (1991). 

10.  Y.  Ovchinnikov  et  al.,  J.  Phys.  B  24*  L539  (1991). 

11.  A.  Sidorov  et  al.,  J.  Phys.  B  24*  3733  (1991). 

12.  K.  Moler  et  al.,  Phys.  Rev.  A45.  342  (1992). 

13.  R.  Grimm  et  al.,  Phys.  Rev.  Lett.  ££,  1415  (1990). 

14.  P.  Hemmer  et  al. ,  Phys.  Rev.  Lett.  ££,  3148  (1992). 

15.  M.  Kasevich  and  S.  Chu,  Phys.  Rev.  Lett.  £2,  1741  (1992). 

16.  C.  Cohen-Tannoud ji  and  W.  Phillips,  Phys.  Tod.  43.*  33  (1990). 

17.  R.  Gupta,  S.  Padua,  C.  Xie,  H.  Batelaan,  and  H.  Hetcalf, 
submitted  to  Opt.  Lett. ,  June,  1993. 

18.  S.  Padua  et  al.,  Phys.  Rev.  Lett.  2£,  3217  (1993). 

19.  We  thank  New  Focus  Corp.  for  the  loan  of  an  EOM. 

20.  S-Q.  Shang  et  al.,  Phys.  Rev.  Lett.  £7,  1094  (1991). 

21.  S-Q.  Shang  at  al.,  Phys.  Rev.  Lett.  ££,  317  (1990). 

22.  G.  Nienhuis  et  al.,  Phys.  Rev.  A44.  462  (1991). 

23.  P.  van  der  Straten  et  al.,  Phys.  Rev.  A47.  4160  (1993). 


9 


( 


•  I  • 


•  • 


t 


4 


FIGURE  CAPTIONS 

Figure  1.  (a)  shows  the  spatial  dependence  of  the  unequal  light 

shifts  (caused  by  either  different  transition  strengths  or 
different  light  intensities)  near  the  two  frequencies  w*  and 
W2 .  (b)  shows  the  atomic  levels  and  definitions  of  the  6's. 

Figure  2.  (a)  shows  the  measured  atomic  beam  profile  for  *  »  0 

and  it.  Here  Si  *  $2  “  46  MHz  (6  *  0)  and  the  intensity  Ii  * 
3.2  mW/cm2  and  I2  *  1.6  mW/cm2  («  Isat  ■  hc/X3r).  A  change 
of  1.0  on  the  vertical  axis  corresponds  to  a  100%  change  of 
the  signal.  (b)  shows  how  the  measured  central  peak  height 
at  large  6  varies  with  #.  The  solid  line  is  a  fit  of  a  sine 
wave  to  the  data. 

Figure  3.  (a)  shows  the  measured  atomic  beam  profile  for  various 

values  of  S  that  produce  VSR  at  velocities  ±6/ 2k.  A  change 
of  1.0  on  the  vertical  axis  corresponds  to  a  100%  change  of 
the  signal.  (b)  shows  the  measured  resonant  velocities  the 

peaks  vs  6.  The  straight  line  is  a  plot  of  v  =  6 /2k.  The 

laser  intensity  is  (0.6)Isat  *  1  mW/cm2  for  each  beam  and  *1 
s  10  MHz. 

Figure  4.  The  spatial  dependence  of  the  light  shifts  for  61  and 
62  with  opposite  signs  and  *  different  from  0  or  n.  The 
arrows  show  the  position  and  the  direction  of  optical  pumping 
between  states  |l>  and  |2>.  Atoms  travel  along  the  thickened 
lines  to  the  right,  and  for  the  case  shown,  the  average  force 
is  clearly  to  the  right.  For  motion  to  the  left,  the  force 
is  still  to  the  right. 

Figure  5.  (a)  and  (b)  show  the  measured  spatial  profile  of  the 
deflected  atomic  beam  for  two  different  signs  of  the  s.  A 
change  of  1.0  on  the  vertical  axis  corresponds  to  a  100% 
change  of  the  signal.  Similar  results  are  obtained  by  moving 
the  retr ©reflecting  mirror.  The  limit  of  velocity  capture 
range  is  evident  in  the  data.  (c)  and  (d)  show  the 
calculated  force  curves.  The  laser  intensity  is  a  1.6  mW/cm2 
for  all  parts,  6*  *  +6  MHz  and  62  *  *3  KHz  for  (a)  and  (c) , 
and  61  «  -3  MHz  and  62  •  +6  MHz  for  (b)  and  (d) . 
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Abstract 
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We  describe  calculations  that  predict  and  analyze  distinctive  quantum 
features  in  the  velocity  distribution  in  ID  laser  cooling  for  light  shift  poten¬ 
tial  well  depths  only  a  few  times  the  recoil  energy.  These  features  arise  from 
transitions  between  energy  bands,  or  even  between  Bloch  states,  in  the  peri¬ 
odic  potential.  They  occur  with  a  <r+  standing  wave,  both  with  and  without 
a  small  B  field,  and  for  lin  X  lin  laser  coding.  We  have  observed  these  fea¬ 
tures  experimentally  in  a  beam  of  metastable  helium  atoms  coded  on  the 

2sS]  *■*  23Pj  transition,  with  velocity  resolution  ~  0.3  recoil. 

PACS  numbers:  32.80.Pj,  42.50.Vk 
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Recent  progress  in  laser  cooling  has  produced  atoms  with  mechanical  energies  less  than 
the  light  shifts  induced  by  the  laser  light  [1,2].  This  means  atoms  may  be  confined  in  the 
periodic  array  of  potential  wells  between  the  nodes  (or  antinodes)  of  optical  standing  waves. 
At  such  low  energies,  the  deBroglie  wavelength  of  the  atoms  is  comparable  to  the  optical 
wavelength.  Motion  in  these  wells  is  thus  quantized:  the  eigenstates  of  atomic  motion 
show  a  band  structure  analogous  to  that  of  electrons  moving  in  the  periodic  potential  of 
a  crystal.  Since  a  theory  of  this  band  structure  in  laser  cooling  was  developed  [3],  there 
has  been  a  keen  interest  in  studying  its  effects  experimentally.  Recently,  the  frequency 
intervals  between  energy  bands  have  been  investigated  with  a  probe  laser  using  four- wave 
mixing  [4],  through  sidebands  in  the  atomic  fluorescence  [5],  and  through  rf  spectroscopy 
[6].  Up  to  now,  however,  there  have  been  no  observations  of  quantum  effects  in  the  velocity 
distribution,  which  is  arguably  closer  to  the  focus  of  laser  cooling. 

Here  we  report  measurements  and  calculations  that,  for  the  first  time,  do  reveal  quan¬ 
tum  effects  in  the  velocity  distribution  of  cooled  atoms.  We  find  not  only  effects  of  tran¬ 
sitions  between  the  bands,  but  also  effects  of  transfer  of  population  between  Bloch  states 
within  these  bands.  Effectively,  we  are  probing  the  velocity  substructure  of  bands  in  the 
periodic  light  shift  potential.  Observation  of  these  features  requires  that  the  velocity  res¬ 
olution  be  smaller  than  the  recoil  velocity  10).  Furthermore,  the  laser-atom  interaction 
rime  must  be  at  least  a  few  times  longer  than  the  optical  pumping  time,  and  the  light 
shift  potential  wells  must  be  shallow  enough  for  only  a  few  bands  to  occur  below  the  tops 
of  the  potential  hills.  These  last  two  conditions  are  met  by  dunce  of  laser  intensity  and 
detuning,  while  the  first  one  is  facilitated  by  cooling  metastable  He  2*$  atoms  (He*)  on 
the  2*S> «-»  23Pa  transition  .  He*  has  a  large  recoil  velocity  Vj»  s  tik/M  -  9.2  cm/s  due 


to  its  small  mass  M  (A  =  2n /k  =  1.083  pm  is  the  laser  wavelength). 

We  discuss  three  cooling  configurations  where  we  have  observed  energy  band  or  Bloch 
state  transition  effects.  The  first  is  a  pure  o+  laser  standing  wave  driving  the  He*  J  = 
1  — »  2  transition,  which  produces  a  dip  in  the  velocity  distribution  near  V  =  0.  The 
dip  is  caused  by  differences  in  the  zero-point  velocity  distributions  in  light  shift  potential 
wells  of  different  depths.  The  second  configuration  is  Magnetically  Induced  Laser  Cooling 
(MILC)  where  a  a+  standing  wave  and  a  magnetic  field  produce  sub-Doppler  cooling  [7]. 
Transitions  between  Bloch  states  within  the  lowest  energy  band  in  the  standing  wave  cause 
a  peak  in  the  velocity  distribution  near  V  =  0  whose  width  is  about  equal  to  Vr.  The  third 
case  occurs  in  polarization  gradient  laser  cooling  with  counter-propagating  laser  beams  of 
orthogonal  linear  polarizations  (lin  ±  lin)  [8].  Here  we  find  dispersion-like  features  in  the 
velocity  distribution  near  V  —  ±V«,  caused  by  net  transfer  of  population  from  states  just 
above  the  lowest  band  gap  to  states  just  below  it. 

The  experimental  apparatus  used  in  our  studies  is  shown  schematically  in  Fig.  1.  He* 
atoms  come  from  a  discharge-excited  nozzle  source,  cooled  by  liquid  JV2.  Through  time-of- 
flight  measurements  we  established  that  the  supersonic  longitudinal  velocity  distribution 
has  an  average  velocity  v  =  1400  m/s  and  a  velocity  spread  A vrus  =  240  m/s  at  our 
standard  operating  conditions.  The  on-axis  source  output  is  —  5  x  10w  He*  atoms  s~* 
sr~l,  while  the  ratio  of  singlet  to  triplet  metastabks  in  the  beam  is  <  0.02.  About  1  cm 
downstream  of  the  nozzle  there  is  a  conical  skimmer,  followed  25  cm  further  by  a  30  pm 
wide  x  7  mm  high  slit  that  reduces  the  flux  to  2  x  10*  He*  atoms  s"1. 

Directly  after  this  slit,  the  atoms  interact  with  two  counter-propagating  laser  beams 
of  equal  intensity.  The  beams  have  a  nearly  Gaussian  intensity  profile  with  a  waist  radius 


ft 


ft 


(*> 


(*) 


I 


> 


ft 


ft  • 


I 


I 


> 


ft 


•  •  •  • 


3 


ft 


w  sb  16  mm  (e-3  intensity  point)  and  tue  apertured  to  a  diameter  of  32  mm.  They  originate 

from  a  home-built,  laser  diode-pumped  CW  LNA-laser  [9].  The  laser’s  frequency  is  locked 

to  a  cavity  that  in  turn  is  locked  to  a  Zeeman-tuned  [10],  weak  RF  discharge  in  helium. 

The  interaction  region  is  surrounded  by  three  orthogonal  pairs  of  Helmholtz  coils  to  create 

a  well-defined  magnetic  field  B  .  The  Earth’s  magnetic  field  was  cancelled  to  ±5  mG  using 

the  Mechanical  Hanle  Effect  on  the  23Sj  — >  23Pi  transition  [11]. 

The  transverse  velocity  distribution  of  the  He*  beam,  modified  by  the  interaction  with 

the  laser  light,  is  measured  with  a  detector  1.9  m  downstream  of  the  interaction  region. 

This  detector  consists  of  a  stainless  steel  plate  placed  behind  a  30  pm  wide  slit  in  a 

moveable,  metal  cage.  The  two  30  pm  slits  separated  by  1.9  m  provide  velocity  resolution  of 

2.7  cm/s  or  0.3  VR.  The  internal  energy  of  the  He*  atoms  can  free  an  electron  from  the  plate 

on  impact,  and  the  resulting  current  (~  10  fA)  is  amplified,  collected  and  recorded.  The 

detector’s  transverse  position  in  the  atomic  beam  is  scanned  with  a  computer  controlled 

* 

stepper  motor. 

We  have  used  two  types  of  density  matrix  quantum  cooling  calculations  [3,12-14] 
to  interpret  the  measurements.  The  density  matrix  p  obeys  the  Liouville  equation 
p  =  ~(i/h)[2f,p]  +  pM,  where  H  includes  internal  energies,  kinetic  energies,  atom-laser 

v 

interactions,  and  magnetic  field  interactions  if  present,  and  pM  expresses  the  spontaneous 
radiative  decay  and  repopulation  of  ground  state  sublevels.  The  calculations  that  best 
model  the  experiment*  use  basis  states  that  are  products  of  internal  atomic  states  and 
free  particle  momentum  eigenfunctions  [3,131;  $  =  |J,m,p  >.  Excited  states  are  explic¬ 
itly  included  so  Doppler  cooling  effects  occur  in  the  results,  and  one  can  also  include  the 

spatially  varying  User  intensity  profile.  These  calculations  match  the  measured  velocity 
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distributions  very  well,  but  do  not  offer  a  clear  picture  of  the  cooling  process. 

More  insight  is  provided  by  our  second  method  because  it  explicitly  uses  eigenstates 
of  the  periodic  light  shift  potential  obtained  by  adiabatically  eliminating  the  excited  state 
[12,14,15].  It  simulates  the  experiment  less  accurately  because  there  is  no  Doppler  cooling 
and  the  laser  profile  is  assumed  constant  Nevertheless,  it  shows  the  key  features  found 
both  in  the  experiments  and  in  the  free  particle  calculations.  The  effective  Hamiltonian 
for  the  ground  state  is  J5fe//  “  (2 S6/L)VgtVtg,  where  $  =  2Q3/rJ(=  1  for  intensity  I  = 
irhc/3A3r),  Q  is  the  single  beam  Eabi  frequency,  T  =  1/t  is  the  decay  rate,  6  =  uWr-watom 
is  the  laser  detuning  from  the  atomic  frequency,  and  1  =  1  +  462/T2.  The  matrix 
includes  Clebsch-Gordan  coefficients  for  the  atom-laser  interaction  (|V^iejl  <  1).  He}j  is 
diagonal  in  the  magnetic  quantum  number  m  for  a*  light  or  for  lin  X  lin  cooling  on  a 
J9  =  5  transition.  The  well  depth  Uo  =  fS6/L ,  where  /  =  2  for  a  o+  standing  wave,  and 
/  =  5/6  for  J  =  1  -*  2  for  lin  X  lin. 

To  help  understand  the  following  effects,  some  background  on  Bloch  states  is  useful. 
When  the  light  shift  potential  is  diagonal,  the  eigenfunctions  satisfy  a  Mathieu  equation. 
Solutions  may  be  written  as  an  expansion  in  momentum  eigenstates  [15,16,3]  = 

Ej  a»J,*exp[t(v+2ji)i:x]  where  n  is  the  band  number  and  v  is  the  Bloch  index:  -1  <  v  <  1. 
Thus  4*(x)  obeys  Floquet's  theorem.  In  Table  I  we  present  the  low  order  wavefunctions 
for  Htjj  —  f/0sina(kx),  when  Ua  is  less  than  the  recoil  energy  Er  =  MVr2/2  =  h2k2/2M. 
In  general,  'i(x)  represents  a  travelling  wave,  except  when  v  s=  0  or  ±1,  where  energy 
gaps  occur.  For  the  few  bands  below  the  tops  of  the  potential  hills  when  the  potential  is 
•hallow,  states  just  below  a  gap  are  more  strongly  localized  in  the  deep  part  of  the  wells 
than  states  just- above  it.  Going  from  the  bottom  to  the  top  of  a  band,  the  wave  functions 


go  from  being  localized  away  from  the  deep  part  of  the  wells  to  being  most  localized  in 
the  deep  part  of  the  wells.  In  a  shallow  potential,  the  lowest  (n  =  0)  band  is  an  exception 
in  that  its  v  =  0  state  is  only  slightly  localized. 

To  explain  our  observations  in  the  clearest  way,  we  use  a  simple  model. 

We  use  a  basis  of  eigenstates  of  the  periodic  light  shift  potential  and  neglect  off-diagonal 
density  matrix  elements  [12]  (except  for  MILC,  where  they  are  always  needed),  and  assume 
that  the  atomic  recoil  from  spontaneous  emission  of  a  *  (a)  photon  is  always  directed 
perpendicular  (parallel)  to  the  laser  axis  (i-axis).  In  a*  excitation  from  m  =  — |  Bloch 
states,  followed  by  it  decay  to  m  =  ]  states,  the  net  change  in  the  Bloch  index  is  then 
Av  as  ±1.  Calculations  based  on  this  simplified  model  preserve  the  features  found  in  more 
exact  calculations. 

We  first  discuss  the  quantum  feature  in  a  a*  standing  wave  (no  polarization  gradients 
or  magnetic  field).  The  data  in  Fig.  2a  show  an  overall  bell  shape  that  arises  from  Doppler 
cooling,  while  the  dip  near  V  =  0  arises  from  purely  quantum  effects.  Fig.  2a  also  shows 
computational  results  from  the  free  particle  basis  (with  Gaussian  laser  spatial  profile),  and 
from  the  periodic  potential  basis  for  a  J  s  |  — *  |  transition. 

The  dip  at  V  =  0  is  a  quantum  feature  that  can  be  explained  even  with  the  sample 
model  Fig.  3a  shows  the  excitation  and  decay  branches,  and  Fig.  3b  shows  the  two 
m  «  potential  wells.  Because  atoms  in  the  shallower  well  have  a  smaller  zero  point 
vibrational  energy  tW  those  in  the  deeper  well,  the  velocity  distribution  in  the  lowest 
band  of  the  shallower  well  is  narrower.  Am  =  1  excitation  followed  by  Am  *  0  decay 
transfer*  atoms  from  the  lowest  band  in  the  shallower  me*]  potential  well  to  the  bwest 
band  in  the  deeper  m  =  +J  well.  This  process  replaces  the  narrow  velocity  distribution 
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in  the  shallow  well  with  a  broader  one.  Fig.  3c  shows  calculated  velocity  distributions 
P(m,n,V)  =  P(-|,0,V)  and  P(f,0,V)  (n  =  0,  1,  ...  is  thv.  band  number),  and  their 
difference  A P{V)  (enlarged  by  a  factor  oi  4),  for  typical  laser  parameters.  AP(V)  has 
a  dip  at  i  =  0  and  two  peaks  at  V  ~  ±1.3V/?.  Thus  this  quantum  effect  is  caused  by 
differences  in  the  velocity  distributions  of  the  lowest  states  of  the  potential  wells  [17].  This 
same  description  holds  for  the  J  =  1  -*  2  transition  in  He*,  and  causes  the  dip  shown  in 
Fig.  2a. 

This  non-Gaussian  population  distribution  is  relatively  short  lived.  For  the  J  =  1  — ♦  2 
transition  in  He*  in  a  o+  polarized  standing  wave,  the  populations  of  the  m  =  -1  and  0 
states  are  depleted  at  the  optical  pumping  rate  =  2ST/L.  For  the  conditions  of  Fig. 
2a,  the  dip  appears  after  ~  1/rs,  which  is  a  few  times  1  /Fp  ~  0.3/is.  However,  Doppler 
cooling  of  m  =  1  atoms  continues  at  the  rate  T o0pp  =  laSSEn/tiTL2.  Since  our  interaction 
time  is  only  ~  20  fi s,  we  do  not  see  the  dip  vanish,  but  our  calculations  indicate  that  it 
disappears  after  about  100  /rs,  which  is  a  few  times  1/Td0PP  ~  20ps. 

The  next  quantum  feature  we  discuss  occurs  when  a  transverse  magnetic  field  is  added 
to  the  <7+  standing  wave  (MILC  [7]).  At  low  intensities  the  two  peaks  discussed  above 
appear  at  approximately  ~  V  =  ±1.3Vr,  but  there  is  also  a  central  peak  of  width  —  V/? , 
as  shown  in  Fig.  2b.  The  population  distributions  in  the  lowest  bands  are  shown  in  Fig. 
4a  vs.  energy.  Fig.  4a  shows  that  the  population  decreases  sharply  with  energy  in  the 
lowest  m  =  -1  band.  No  central  peak  occurs  when  the  populations  of  the  Bloch  states 
in  this  band  are  uniform,  as  they  nearly  are  for  the  pure  a*  case.  We  conclude  that  the 
variation  in  the  population  over  the  band  is  responsible  for  the  peak  in  P(V)  at  V  =  0  in 
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Fig.  2b. 
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As  discussed  in  Ref.  (13],  MILC  is  a  cyclic  process  of  Am  =  +1  optical  pumping 
and  Am  =  ±1  magnetic  field  mixing.  (Hence  in  MILC,  the  peaks  at  V  =  ±1.3V*  are 
continually  regenerated.)  In  Fig.  4b,  transitions  between  Bloch  states  (labeled  A  through 
D)  in  the  lowest  bands  (n  =  0)  of  the  m  =  i|  potentials  are  shown  by  straight  and 
curved  (single)  arrows,  respectively.  The  numbers  (0.87,  etc.)  give  the  relative  pumping 
transition  rates  between  states  |m, n, v  >  as  Tad  ~  I  <  — 1,0, 0|stn(A:e)||,0, 1  >  |2  =  0.80 
and  Tbc  —  I  <  -|,0,l|stn(ir)||,0,0  >  |J  =  0.87.  In  the  deeper,  m  =*  +§  well,  all  states 
are  comparably  localized  near  the  high  intensity  region,  but  in  the  shallow  m  =  -|  well, 
the  n  =  0  wave  functions  become  more  strongly  localized  in  the  high-intensity  part  of  the 
well  as  v  goes  from  0  to  1.  Therefore  Tbc  is  greater  than  Tad-  The  pumping  rate  from 
m  =  -J  to  m  w  i  for  v  -  0  (point  A)  to  v  -  1  (D)  is  weaker  than  the  pumping  rate 
from  i/  -  1  (B)  to  i/  a  0  (C)  because  atoms  at  B  experience  a  higher  average  intensity 
than  those  at  A.  By  contrast,  the  magnetic  field  mixing  rate  between  the  two  v  =  1  states 
(BD)  is  very  nearly  equal  to  that  for  v  =  0  (AC).  The  net  result  of  this  cyclic  process  is 
an  accumulation  of  atoms  in  the  state  with  the  lowest  pumping  rate,  i.e.,  state  A,  (double 
arrow  in  Fig.  4b).  As  seen  from  Table  I  and  Fig.  4c,  atoms  with  v  m  0  have  a  large 
probability  to  be  at  V  =  0.  Thus,  a  central  peak  occurs  because  of  a  two-step  transfer 
of  atoms  within  the  lowest  m  =  -|  energy  band,  and  is  a  purely  quantum  feature.  Fbr 
deeper  potentials  with  more  bands  below  the  tops  of  the  potential  hills,  this  quantum 
feature  evolves  into  the  usual  MILC  cooling  peak. 

The  third  quantum  feature  we  discuss  in  ID  laser  cooling  occurs  in  lin  J.  lin  polarization 
gradient  cooling  at  low  intensities.  For  Uo  **  Er,  we  find  that  the  velocity  distribution 
has  strong  dispersion  shapes  centered  at  V  =  ±V*,  which  can  also  be  understood  in 
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terms  of  transitions  between  Bloch  states  within  the  low-lying  energy  bands.  Calculations 
and  experimental  data  on  this  effect  are  shown  in  Fig.  2c.  To  obtain  the  best  match 
between  experiment  and  theory  in  Fig.  2c,  the  optimum  laser  intensity  parameter  for  the 
calculations  was  found  to  be  S  =  1.2  as  compared  with  the  nominal  experimental  value  of 
2.7.  This  discrepancy  is  unexplained  at  present. 

In  lin  1.  lin  cooling,  each  optical  pumping  step  of  the  cyclic  transfer  from  m  =  —  | 
to  m  as  i  and  back  results  in  some  energy  less.  With  shallow  potentials,  the  lowest 
quantum  states  produce  anomalies  in  the  cooling  process.  We  focus  on  the  optical  pump¬ 
ing/radiative  decay  transition  from  the  m  =  —  |,n  =  1  band  to  the  m  —  |,n  =  0  band. 
Atoms  in  Bloch  states  |m, n,  v  >=  |  -  |,  1,  ;$  1  >  are  pumped  to  ||,0, ~  0  >  by  <r+  pho¬ 
tons,  then  to  |  -  0,1$  1  >  by  a~  photons  (Fig.  5).  FVom  |  -  |,0,<  1  >,  atoms  are 

pumped  to  higher  states  at  rates  smaller  than  for  these  processes,  so  there  is  a  net  loss  of 
j  -  5, 1,  £  1  >  atoms  and  a  net  accumulation  of  |  -  |,  0,  1  >  atoms.  For  |  -  5, 1,  1  >, 
the  most  probable  |V|  is  slightly  >  V*,  while  for  |  —  1  >,  the  most  probable  |V|  is 

slightly  <  V/i,  as  seen  from  the  velocity  distributions  in  Fig.  5b  or  the  (1,0.8)  and  (0,0.8) 
functions  in  Table  I.  This  two-step  transition  accounts  for  the  dip  at  V  >  V/j  and  the  peak 
at  V  ;$  V/j,  and  reflects  the  differing  character  of  Bloch  states  just  above  and  just  below 
the  lowest  band  gap. 

In  summary,  the  experiments  and  calculations  reported  here  show  that  in  order  to 
understand  laser  cooling  on  the  level  of  one  recoil  velocity,  one  must  understand  the 
redistribution  of  population  Bloch  states  in  the  periodic  light  shift  potential. 

This  work  was  supported  by  NSF,  ONR,  AFOSR,  the  US  Dept  of  Education,  and  by 
a  grant  of  computer  time  from  the  Cornell  National  Supercomputer  Facility. 
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TABLES 

TABLE  I.  Low  order  expressions  for  some  unnormalized  Mathieu  functions,  labeled  by  (n,  u) 
(adapted  from  Ref.  [16]).  Here  u  =  Uo/Er  and  z  =  kz.  Energies  (third  column)  are  relative  to 
the  potential  mean. 


Ob*') 

* 

E/Er 

0,0 

1  -  (u/8)  cos(2z) ... 

— ua/32  ... 

0,0.8 

exp(0.8u)  -  (5u/16)exp(-1.2iz) ... 

0.64  -  0.09ua  ... 

0,1 

sin(r)  -  (ti/32)sin(3r)  ... 

1  -  u/4  ... 

1,1 

cos(z)  -  (u/32)cos(3r)  ... 

1  +  u/4  ... 

1,0.8 

exp(-1.2iz)  +  (5u/16)exp(0.8tz) ... 

1.44  +  0.07ua  ... 

FIGURES 


FIG.  1.  Schematic  of  the  experimental  apparatus. 

FIG.  2.  Experimental  data  (dots)  and  free  particle  basis  computational  results  (solid  lines)  for 
the  velocity  distribution  of  He*  atoms  with  (a)  a  a*  standing  wave ,  Sma*  =  1.65, 6  =  — 2r;  (b)  a  <r+ 
standing  wave  plus  50  mG  B  field,  S^tr  -  4.5, 6  =  -4F;  (c)  lin  J.  lin  cooling,  5mat  —  2.7,5  =  -8F. 
For  (c),  calculations  were  performed  with  Snax  =  1.2.  The  laser  intensity  has  a  Gaussian  spatial 
distribution.  Dashed  lines  give  results  of  periodic  potential  basis  calculations  for  the  simpler  case 
of  J  =  1/2  -*  3/2,  with  adjusted  uniform  laser  intensity. 

FIG.  3.  Analysis  of  the  quantum  feature  in  Fig.  2a,  using  a  J  =  1/2  — *■  3/2  transition,  (a) 
Excitation  and  decay  scheme,  with  Clebsch-Gordan  factors  and  m  values  as  shown;  (b)  One  period 
of  the  light  shift  potential  for  m  =  ±1/2,  with  the  lowest  energy  bands  indicated;  (c)  Velocity 
distributions  for  the  lowest  energy  bands,  and  their  difference. 

FIG.  4.  Analysis  of  the  quantum  feature  in  Fig.  2b  (M1LC).  (a)  Calculated  populations  in  Bloch 
states  (discretized  to  a  mesh  of  —  0.1)  for  the  lowest  energy  bands  for  m  =  ±1/2.  (b)  Low-lying 
energy  bands  tor  m  =  ±1/2,  with  optical  pumping  (straight  arrows)  and  magnetic  field  coupling 
(curved  single  arrows)  transition  elements,  v  is  the  Bloch  index.  The  double  arrow  shows  the  net 
population  transfer,  (c)  Velocity  distribution  functions,  P(V),  for  the  m  =  -1/2, n  =  0,i>  =  0  and 
1  Bloch  states,  initial  and  final  states  of  the  two-step  transfer  process. 

FIG.  5.  Analysis  of  the  quantum  feature  in  lin  1  lin  cooling  in  Fig  2c,  in  terms  of  a 
J  b  1/2  -♦  3/2  transition,  (a)  Low  lying  energy  bands  for  m  =  ±1/2,  with  arrow  showing 
important  transitions  between  Bloch  states,  (b)  P(V)  functions  for  the  initial  and  final  Bloch 
states  in  the  two-step  transfer  shown  in  (a). 
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Abstract 

We  have  demonstrated  slowing  of  a  rubidium  atomic  beam  by  isotropic 
monochromatic  light.  The  results  agree  with  a  model  calculation,  thus  al¬ 
lowing  its  use  for  designing  an  isotropic  light  slower.  The  large  hyperhne 
splittings  of  rubidium  lead  to  natural  multi-frequency  slowing,  which  is  also 
included  in  our  model. 

PACSnr.  32.80.Pj,  42.50.Vk 

I.  INTRODUCTION 

The  ultimate  goal  for  atomic  beam  slowing  and  coding  is  the  compression 
of  the  entire  velocity  distribution  of  atoms  in  a  beam  as  they  emanate  from 
a  source  to  one  narrow  peak  at  a  desired  low  velocity.  Although  this  goal  is 
far  from  being  reached  at  present,  atomic  beam  slower*  are  commonly  wed 
for  atomic  docks  and  for  loading  of  magneto-optical  traps  and  optical  mo¬ 
lasses.  Improvement*  to  atomic  baam  slower*  are  therefore  of  current  interest. 
Recently  slewing  and  cooling  by  white  light  {1}  and  isotropic  monochromatic 
light  (2]  have  been  demonstrated  to  be  promising  new  techniques. 

In  isotropic  User  light,  the  rasonsnee  frequency  u,/  of  an  atom  moving  at 
a  velodiy  v  will  be  matched  by  the  Doppler-shifted  User  frequency  when  the 
wavevector  I  of  the  light  make*  an  nng)e  #  with  the  ntomic  velocity  $; 


H/  * t  mi  -  keens#, 


where  vi  it  the  laser  frequency  and  t  (/)  indicates  the  initial  (final)  state.  For 
red  detuned  light  (uj  -  w,-/  <  0)  this  resonance  condition  requires  9  >  x/2, 
meaning  that  the  recoil  due  to  the  absorbed  photons  opposes  the  atomic  mo¬ 
tion  and  slows  the  atoms.  Subsequent  spontaneous  emission  does  not  exert 
a  force  on  the  atoms  on  average.  As  the  atoms  are  decelerated,  they  absorb 
light  from  a  increasing  angle  6  (angle  tuned),  until  the  maximum  value  of 
6  (=  »)  is  reached.  Because  the  light  is  isotropic,  the  atomic  motion  it. 
rectly  opposed  irrespective  of  its  direction,  while  for  Zeeman  (magnetically 
tuned)  [3]  and  chirped  (frequency  tuned)  [4]  slowing,  only  the  longitudinal 
velocity  component  of  the  atomic  motion  is  opposed.  This  major  advan¬ 
tage  prevents  the  atomic  beam  from  being  transversely  expanded  (apart  from 
spontaneous  emissions),  in  contrast  to  Zeeman  tuned  and  chirped  slowing. 

In  this  article,  we  describe  the  application  of  the  isotropic  slowing  tech¬ 
nique  to  a  rubidium  beam  in  order  to  test  experimental  techniques  and  a 
numerical  model  for  the  slowing  process.  The  resulting  slowed  beam  of  cold 
atoms  is  enhanced  by  natural  multi-frequency  slowing  [2],  for  which  the  ru¬ 
bidium  atom  has  an  appropriate  hyperfine  splitting  of  the  excited  state.  The 
carefully  tested  model  is  then  used  to  calculate  isotropic  slowing  over  an  ex¬ 
tended  length. 

11.  EXPERIMENT 

A  schematic  view  of  oar  apparatus  is  shown  in  Fig.  1.  A  thermal  beam 
of  Rb  is  collimated  to  —10  mrad  by  an  aperture  (Al)  suspended  from  a  cold 
trap,  and  travels  through  a  slower  made  from  the  diffusely  reflective  material 
Spectraloa  ($].  The  20  cm  long  cylindrical  slower  has  an  inner  diameter  of  3 
mm  and  outer  diameter  o£  2  cm.  The  reflectivity  r  of  the  material  is  specified 


by  the  manufacturer  to  be  99.1%.  A  beam  from  a  Ti:sapphire  laser,  split  into 
two  branches  by  a  beam  splitter  (BS),  enters  the  tube  via  two  1  mm  diameter 
holes  situated  5  cm  from  each  cylinder  end.  The  laser  light  has  a  line  width  of 
<  1  MHz  and  is  side-locked  to  the  Doppler  broadened  bSi/2  -»5P3/2  transition 
of  Rb  near  A  =  780  nm  with  an  accuracy  of  ±  50  MHz  using  absorption  in  a 
vapor  cell.  Saturated  absorption  is  used  to  calibrate  the  Doppler  broadened 
absorption  spectrum  frequency  scale.  An  electro-optical  modulator  (EOM)  [6] 
divides  the  laser  power  between  the  main  laser  frequency  wj  =  2r  x  (Fig. 
2)  and  two  side  bands  at  W|  ±  2v  x  2.91  GH2.  One  sideband  (1/3)  counteracts 
optical  pumping  into  a  single  ground  state  hyperfine  sublevel,  while  the  othei 
is  not  used.  Alternatively,  we  use  the  Titsapphire  laser  and  a  Sharp  model 
LT025  diode  laser  (DL1)  to  provide  two  frequencies,  instead  of  the  EOM. 

The  laser  beam  from  another  Sharp  model  LT024  diode  la-er  (DL3)  is 
directed  ai  an  angle*  of  30"  with  the  atomic  beam  to  the  determine  the  longi¬ 
tudinal  velocity  distribution  by  measuring  the  intensity  of  the  emitted  fluo¬ 
rescence  as  a  function  of  the  laser  frequency.  The  saturated  absorption  signal 
from  this  laser  also  serves  as  a  calibration  of  its  frequency.  A  beam  from  a  fre¬ 
quency  locked  laser  (DL2)  perpendicular  to  the  atomic  beam  cancels  velocity 
selective  optical  pumping  effects. 

Two  concave  spherical  mirrors  (CM)  collect  the  fluorescence  (7,8)  and  fo¬ 
cus  it  on  a  hybrid  photodiode-amplifier  detector  [9]  (FD).  Stray  laser  light 
emanating  from  the  Spectral  an  tube  U  shielded  from  the  detector  (aperture 
A2).  We  chop  the  tide  bands  of  the  EOM  at  500  Hz  and  subtract  the  gated 
signal  from  the  gated  background  to  reduce  the  noise  markedly. 
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III.  A  MODEL  CALCULATION 


We  assume  that  the  initial  longitudinal  velocity  distribution  of  the  atomic 
beam  is  thermal  and  given  by  [10] 

/(v)  oc  t>3e"wvS^J*®T, 


where  kg  is  the  Boltzman  constant  and  T  is  the  oven  temperature.  The  force 
on  the  atoms  can  be  written  as 


f(*)  =  J  —  DeJ2«jCejp)(iil, 


where  the  subscript  g  (e)  indicates  the  ground  (excited)  hyperfine  levels,  H 
indicates  the  steady  state  population  and  0  indicates  the  direction 

from  which  an  atom  absorbs  light.  The  angle  0  is  defined  relative  to  the 
direction  of  the  atomic  motion.  The  projection  of  the  recoil  momentum  of 
the  atom  p  along  its  direction  of  motion  is  given  by 


p  s  hkcosB 


The  transition  rate  [11]  from  an  initial  (i)  state  to  a  final  (f)  state  is  given  by 


c./  =  (2A  +  l){2f)  +  lH 


h  1  Fi 
Ji  /  Ji 


1 

**"  * 
r 


Here,  the  initial  (final)  state  can  be  either  the  ground  or  excited  state.  The 
lifetime  r  of  the  excited  Jt  state  (for  the  5*P yi  state  of  wRb)  is  27  ns  (12], 
/  *  5/2,  Jt  «  1/2,  J,  a  3/2, =  2,3, F,  a  1, 2,3,4  and  {  }  are  the  usual 
6-j  symbols  (13].  The  photon  scattering  rate  &/(*,t>)  is  proportional  to  the 
convolution  of  the  User  frequency  distribution  /  inside  the  tube  and  the  atomic 
line  shape  S,/(u>)  associated  with  the  transition  from  the  initial  to  the  final 
state  [14] 


■=  J  £  n. 


I 


where  T  (=  2*xt>  MHz)  is  the  atomic  line  width  for  85 Rb  and  p(z)i:  the  light 
energy  density  at  position  x.  In  the  atomic  line  shape 

S'jiu>)  =  (u  -J}?  +  r» 

Uij  =uj-u,  denotes  the  frequency  difference  between  the  hyperfine  levels  F / 
and  Fj.  Here,  I  is  the  normalized  laser  frequency  distribution  in  the  atomic 
rest  frame 

j  Pj 

/«/(».«.*)  =  ^2  +  k„CM$y  +  ff*1* 

where  fj  is  the  laser  line  width.  The  energy  density  p(z)  is  derived  by  an 
iterative  procedure  [15]. 

For  a  Lambertian  reflector,  which  obeys 

J(a)  =  J0coxQ, 

where  Jo  is  the  incident  laser  power  and  J{a)  the  reflected  power  at  an  angle  a 
with  respect  to  the  normal  of  the  reflector  surface,  a  straightforward  derivation 
for  a  cylindrical  geometry  yields 

for  the  power  distribution  Pm  of  the  n-th  reflection  of  the  laser  light  in  the 
cylinder.  The  initial  power  distribution  is  given  by  Pq.  For  example,  with  one 
entrance  hole  at  x  =  0  the  total  laser  power  P  =  Po(0)  and  Pq(x  jf  0)  «  0. 
The  length  of  the  cylinder  is  l  and  the  reflectivity  r.  We  have  assumed  that 
the  power  distribution  is  symmetric  in  4-  The  total  power  distribution  may 
be  obtained  by  summing  over  nil  reflections,  P(x)  »  /*(*).  We  can 

immediately  calculate  the  energy  density  bom  [10], 

p(*)  =  4P(*)/e 

by  assuming  that  the  power  distribution  is  isotropic.  We  have  verified  ana¬ 
lytically  that  for  an  infinite  cylinder,  for  which  the  surface  is  a  homogeneous, 


# 
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Lambertian  source  the  power  distribution  is  isotropic.  These  approximations 
are  motivated  by  the  small  Tadius  over  length  ratio  for  the  cylinder. 

In  the  force  expression  the  steady  state  populations  He,H«  of  the  hyperfine 
levels  as  a  function  of  the  position  z  are  calculated  by  substituting  Rij(z,v) 
in  the  rate  equations  [11,16).  Obtaining  from  the  calculated  force  the  final 
velocity  of  the  atom  vj  for  each  initial  velocity  v,,  we  obtain  the  final  velocity 


distribution, 


expressing  the  conservation  of  atoms. 


IV.  RESULTS 


To  test  the  effective  laser  power  in  the  tube  we  measured  and  calculated 
the  fraction  of  atoms  optically  pumped  from  the  F,=3  into  the  Fa~2  state 
as  a  function  of  laser  power  injected  into  the  Spectralon  tube.  We  chose  this 
fraction  to  be  the  ratio  of  the  maxima  of  the  atomic  velocity  distribution  with 
and  without  laser  power,  and  plotted  the  result  in  Fig.  3.  The  calculation  is 
the  result  of  the  time  dependent  rate  equations  [16]  with  the  isotropic  photon 
scattering  rate  (Eq.  6)  for  r&99%.  For  an  oven  temperature  of  T=  130°C 
the  density  of  rubidium  results  in  an  estimated  mean  free  path  of  the  photon 
of  6  cm,  while  for  an  oven  temperature  of  T=  1?0®C  the  mean  free  path  is 
only  0.6  on.  On  the  other  hand  a  photon  travels  on  average  2.5  cm  before 
being  absorbed  by  a  tube  with  a  92%  reflectivity  (see  below)  and  3  mm  di¬ 
ameter.  Therefore  we  attribute  the  decrease  in  the  optical  pumping  efficiency 
with  increasing  temperature  both  to  the  optical  thickness  of  the  atomic  beam 
and  to  the  strongly  Doppler-shifted  fluorescence  from  all  these  atoms  that 
is  prevented  from  escaping  by  multiple  reflections  in  the  Spectralon.  AH  ve¬ 
locity  distributions  presented  below  were  taken  with  an  oven  temperature  of 


•  • 


T=  130°C  where  the  Atomic  beam  is  sot  optically  thick.  The  optical  pump¬ 
ing  efficiency  test  is  thus  used  as  an  atomic  calibration  of  the  energy  density 
inside  the  Spectralon  tube. 

Still,  roughly  a  factor  of  4  discrepancy  exists  between  the  measured  and 
the  calculated  optical  pumping  efficiency.  An  additional  experimental  test 
revealed  the  source  of  the  problem.  Fig.  4  shows  the  result  of  a  measure¬ 
ment  of  the  relative  light  power  that  has  diffused  to  the  outside  of  the  tube 
as  a  function  of  position  along  it.  The  experimental  result  for  a  Rb-exposed 
Spectralon  tube  before  and  after  cleaning  were  compared  with  the  result  cal¬ 
culated  from  equation  10.  The  results  for  the  relative  light  power  were  scaled 
appropriately.  The  drop  in  reflectivity  of  the  exposed  tube  corresponds  to  a 
factor  of  4  drop  in  intensity.  In  accordance  with  the  atomic  calibration  of  the 
energy  density,  we  used  a  factor  4  to  correct  the  laser  power  in  all  subsequent 
calculations  for  the  velocity  distributions. 

The  measured  fluorescence  spectrum  was  used  to  determine  the  velocity 
distribution,  which  has  contributions  from  each  of  the  three  accessible  e* cited 
hyperflne  states.  A  given  frequency  of  the  detection  laser  corresponds  to 
three  different  detunings,  which  in  their  turn  via  £,/  =  kvtjcos9  correspond  to 
atoms  with  three  different  lonptudinal  velocities.  In  our  case,  the  repumping 
detection  laser  (DL2)  was  strongly  saturating,  while  the  frequency  ramped 
detection  laser  (DL3)  was  operated  in  the  linear  range,  causing  the  relative 
strengths  of  the  different  contribution*  to  be  proportional  to  «,/.  Taking  these 
three  corrections  into  account  in  our  calculated  distributions  gives  excellent 
agreement  with  the  experimental  mult  for  the  initial  velocity  distribution 
(*1*  5). 

In  Fig.  6  the  difference  between  the  final  and  initial  velocity  distribution 
is  shown.  A  negative  value  of  this  difference  indicates  that  atoms  have  been 
removed,  while  a  positive  value  shows  an  increase  in  the  number  of  atoms.  Fig. 
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6a  shews  that  atoms  around  400  m/s  have  been  slowed  to  velocities  around 
150  m/s.  The  experimental  detuning  dj  (Fig.  2)  is  *190  ±  50  MHz  relative  to 
the  Ft= 3  to  Ft=4  transition  and  the  detuning  used  in  the  calculation  is  -210 
MHz.  The  EOM  is  driven  at  2.91  ±  0.01  GHz,  which  means  that  the  detuning 
6^  relative  to  the  Ft-2  to  Fe=3  transition  is  also  >190  MHz.  The  total  laser 
power  P  =  100  mW,  and  the  ratio  0  of  each  EOM  sideband  over  carrier  is 
0  =  1.  The  vertical  scale  is  normalized  to  the  height  of  the  maximum  of  the 
initial  atomic  velocity  distribution,  so  that  the  peak  of  slowed  atoms  (Fig. 
6a)  contains  ~2%  of  all  the  incident  atoms. 

The  acceleration  in  the  velocity  distribution  shown  in  Fig.  6b  has  been 
measured  under  the  6ame  experimental  conditions  as  in  Fig.  6a,  except  for  a 
positive  190  MHz  detuning.  The  presence  of  the  double  dip  structure  can  be 
mostly  attributed  to  the  Ft= 3  to  Ft=4  and  to  the  F$=Z  to  F*= 3  transitions. 
The  Doppler  resonance  condition  (Eq.  1)  is  met  for  the  Fa= 3  to  F,=4  transi¬ 
tion  at  170  m/s  and  for  the  Fa= 3  to  F,s= 3  transition  at  280  m/s,  which  agree 
with  the  positions  of  the  dips. 

For  negative  detuning  the  Doppler  resonance  condition  (Eq.  1)  is  met  lor 
the  Ff -3  to  F,=4  transition  at  170  m/s  and  for  the  F$= 3  to  Ff=3  transition 
at  60  m/s.  Since  the  interaction  time  is  longer  at  low  velocities,  the  slow¬ 
ing  process  is  more  efficient  than  the  acceleration  process,  and  the  two  peaks 
merge  into  one.  This  is  natural  multi-frequency  slowing  and  is  especially  ef¬ 
fective  because  of  the  appropriate  hyperfiae  splitting  of  Rb.  It  can  not  be 
exploited  in  sodium  (2),  where  the  hyperfine  splitting  is  too  small-  Instead 
one  mat*  assist  the  angle-toned  slowing  with  a  stepwise  change  in  laser  fre¬ 
quency  (multi-frequency  slowing).  The  large  excited  state  hyperfine  splitting 
of  rubidium  provides  .his  effect  n  *  -t  r.-Oy. 

The  results  in  Fig.  oc  and  6d  have  been  measured  for  a  ratio  of  EOM  side¬ 
band  over  carrier  0  *  0.02.  For  positive  detuning  the  resonance- condition  (at 
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170  m/s)  can  only  be  met  by  the  cyclic  F$=3  to  Fe= 4  transition.  Therefore 
a  weak  repumping  sideband  is  sufficient  to  maintain  the  acceleration  process. 
However,  for  negative  detuning  the  f^=3  to  Fe=3  (and  Fg=3  to  Ft- 2)  tran¬ 
sition  can  excite  atoms  at  170  m/s  due  to  angle-tuning.  This  means  that  a 
weak  repumping  sideband  is  not  sufficient  to  maintain  the  slowing  process. 

Finally,  having  established  confidence  in  the  numerical  model  by  the  agree¬ 
ment  with  the  measurements,  we  show  the  result  of  calculations  for  a  20, 60 
and  100  cm  long  slo^ers  with  a  500  mW  input  power  (Fig.  7).  Note  that 
the  fraction  of  atoms  s/ound  50  m/s  has  increased  by  a  factor  of  100  over  the 
initial  distribution.  Because  this  strong  increase  shows  promise  for  important 
applications  of  this  technique,  we  are  currently  building  an  extended  length 
slower. 
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Figure  captions 


Fig.  1.  A  schematic  view  of  the  apparatus.  An  explanation  of  its  com¬ 
ponents  is  given  in  the  text. 


Fig.  2.  A  scheme  of  the  relevant  energy  levels  of  “Rb. 


Fig.  3.  Fraction  of  atoms  in  the  Ff=3  state  after  the  atoms  pass  through 
the  slower  for  oven  temperatures  of  170°  (A)  and  130°  (□)  are  given.  The 
lino  is  the  calculated  result  from  the  rate  equations. 


Fig.  4-  Laser  power  measured  as  a  function  of  position  along  the  external 
wall  of  the  Rb-exposed  Spectralon  tube  of  10  cm  length  before  (o)  and  after 
(A)  cleaning.  The  calculated  power  is  given  for  r  =995?  (•*•)*  98%  ( — )>  965? 

( - )  and  925?  ( - — ),  The  results  are  scaled  appropriately.  The  laser 

beam  enters  the  tube  at  the  aero  position. 


Fig-  5.  The  measured  and  calculated  initial  velocity  distribution. 


Fig.  6.  The  difference  between  the  slowed  and  initial  atomic  veloc¬ 
ity  distribution.  The  experimental  and  calculated  results  show  agreement. 
The  total  laser  power  is  P  »  100  mW.  (a)  Natural  multi-frequency  slowing, 


6  =  -130  ±  50  MHz,  0  =  1.  (b)  Acceleration  showing  two  dips,  6  =  190  ±  50 
MHz,  0  »  1.  (c)  Slowing  disappears,  d  =  -190  ±  50  MHz,  0  =  0.02.  (d) 
Acceleration  remains,  d  =  190  i  50  MH2,  0  =  0.02. 

Fig.  7.  Calculated  initial  (solid  line)  and  final  atomic  velocity  distribu¬ 
tions  for  d  =  -75  MHz,  P  =  500  tnW,  r=92%,  0  *  0.5  and  the  EOM  is 

driven  at  2.97  GHz.  The  lengths  of  the  slower  are  20  cm  (•••)»  60  cm  ( - ) 

and  100  cm  ( - )•  The  largest  slowed  peak  at  50  m/s  contains  -25%  of 


all  incident  atoms. 


velocity  (m/s) 
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